UNCLASSIFIED 


AD  NUMBER 


AD908907 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  only;  Test  and  Evaluation;  JAN 
1973.  Other  requests  shall  be  referred  to 
Air  Force  Armament  Laboratory,  DLDL,  Eglin 
AFB,  FL  32542. 


AUTHORITY 


AFATL  ltr,  10  Oct  1975 


THIS  PAGE  IS  UNCLASSIFIED 


THIS  REPORT  HAS  SEEN  DELIMITED 
mV  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOB  DIRECTIVE  5200 . 20  AND 
NO  RESTRICTIONS  ' ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED, 


AFATL-TR-73-15 


MONTE  CARLO  ANALYSIS  OF  S-CURVE 
AND  ROLL-THROUGH-ZERO  BOMBLET 
DISPERSION  CHARACTERISTICS 


ALPHA  RESEARCH,  INC. 


-  a  It. 


dot?  T«* 


TECHNICAL  REPORT  AFATL-TR-73-15 
JANUARY  1973 


D  D  C 


Distribution  limited  to  U.  S.  Government  agencies  only; 
this  report  documents  test  and  evaluation;  distribution 
limitation  applied  January  1973  .  Other  requests  for 
this  document  must  be.  referred  to  the  Air  Force  Armament 
Laboratory  (uLDL)  ,  Eglin  Air  Force  Base,  Florida  32542. 


AIR  FORCE  ARMAMENT  LABORATORY 

AIR  FORCI  SYSTIMS  COMMAND  •  UNITKD  STATU  AIR  FORCl 


EGLIN  AIR  FORCE  BASE,  FLORIDA 


Best 

Available 

Copy 


Monte  Carlo  Analysis  Of  S-Curve  And  Roll-Through-Zero 
Bomblet  Dispersion  Characteristics 


f f 1  i  st  r i but  ion  1 1 n 
report  doer 
j  1  !  iti  1 1  ;it  i  on  app  1  1 1 
j  t !;  i  >  bo*,  merit  rnr 
j  i.aborat  orv  p  lab., 


domes  E.  Brunk 


*  t  o* 

to 

i.  s,  . 

uvern-sent  .n 

one  i  o 

\ - 1 

s  riUV;  1 

•tent 

S  t  C 

it  arul 

ova  Uiat  ion ; 

■  1 1  •-  r  r 

1 

>  1-1  ■ t  J  I  vf* 

-.!  .1 

•I.'iuan 

.  uthi-r  n 

ouost 

s  t‘o  r  • 

•  t  i 

t*  :  c 

!  erred 

to  the  Air  1 

ot\o 

Amir.enf  ; 

.  1 

1  i  n 

A  i  r  1  ( 

roe  ti-i so.  1  I 

orida 

CcjH’  cr?nf!abi«  to  DOC  dr *$  not 
pmirji  fu:iy  lep.-wdccaao 


7  n~~V^iPiD  >-?  P  1 7"rf  " 


FOREWORD 


This  final  report  documents  work  accomplished  during  the  period 
4  February  1971  through  15  November  1972  by  Alpha  Research,  Inc.,  Santa 
Barbara,  California,  under  Contract  No.  F08635-71  -C-0089  with  the  Air 
Force  Armament  Laboratory,  Eglin  Air  Force  Base,  Florida.  Monitoring 
the  program  for  the  Armament  Laboratory  were  Messrs.  M.  Bouffard  and 
2nd  I.t.  Paul  Mayer  (DLDL).  The  principal  investigator  for  the  contractor 
was  Mr.  James  E.  Brunk. 

Digital  computer  service,  flight-test  support  service,  and  model 
design  and  fabrication  service  in  support  of  this  effort  were  provided  by 
the  Air  Force  Armament  Laboratory  and  Armament  Development  and  Test 
Center.  Eglin  Air  Force  Base,  Florida. 

Additional  wind  tunnel  test  support  was  provided  by  the  4T  Projects 
Branch,  PWT,  Arnold  Engineering  Development  Center,  Arnold  Air  Force 
Station,  Tennessee. 

This  technical  report  has  been  reviewed  and  is  approved. 


Director,  Guns  h  Rockets  Division 


ABSTRACT 


The  flight  dynamics  and  dispersion  characteristics  of  candidate 
S-Curve  and  roll -through-zero  bomblets  are  evaluated  by  use  of  the  Monte 
Carlo  method.  Detailed  analyses  of  bomblet-cluster  breakup  and  bomblet 
configurational  and  aerodynamic  asymmetries  are  used  to  establish  the 
statistical  input  data  required  for  the  Monte  Carlo  analyses.  Improved 
aerodynamic  data  packages  for  both  the  S-Curve  and  roll-through- zero 
bomblets  are  provided,  and  results  of  additional  wind  tunnel  tests  are  dis¬ 
cussed.  Monte  Carlo  simulations  show  that  for  realistic  flight  environ¬ 
ments.  both  the  S-Curve  and  roll-through- zero  bomblets  can  achieve  large 
and  uniform  impact  patterns.  For  a  representative  dispenser  HOB  of  2000 
feet, opening  velocity  of  900  feet/second,  and  45  degrees  flight  path  angle, 
impact  pattern  widths  of  800  and  600  feet  are  computed  for  the  S-Curve  and 
roll-through- zero  type  bomblets,  respectively.  The  effects  of  transonic 
and  supersonic  delivery,  large  static  mass  unbalance,  nose-roughness 
asymmetry,  and  intentional  fin  cant  on  the  S-Curve  bomblet  flight  character¬ 
istics  and  dispersion  are  investigated  in  detail.  Flight  test  dispersion  data 
for  the  S-Curve  bomblet,  as  obtained  from  air-gun  launched  models,  is 
compared  with  results  from  analytical  simulations.  A  dynamic  wind  tunnel 
model  support  system,  designed  for  testing  the  S-Curve  bomblet,  is  de- 
scri  bed. 
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CSF, 
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SECTION  I 


INTRODUCTION 

One  of  the  most  effective  means  of  dispersing  a  cluster  bomb  unit  is 
through  the  use  of  aerodynamic  self-dispersing  bornblet  munitions.  With  this 
approach  the  energy  requirements  for  dispersal  are  not  only  derived  from 
the  kinetic  and  potential  energy  of  the  aircraft  and  dispenser,  but  more  sig¬ 
nificantly,  the  dispersal  force  is  applied  over  the  entire  bornblet  flight 
period.  The  advantage  of  aerodynamic  dispersion  increases  with  increasing 
delivery  Mach  number. 

While  self-dispersing  magnus  -  rotor -type  bomblets  have  received 
considerable  attention  and  utilization  in  recent  years,  other  types  of  aero¬ 
dynamic  sel f- di spe r  si  ng  bomblets  have  yet  to  be  introduced  in  quantity.  Two 
new  types  of  sol  f  -  di  spe  r  s  i  ng  bomblets  appear  particularly  promising,  namely, 
the  S-Curve  typo  bornblet  and  the  r  oil -th  rough  -  ze  ro  type  bornblet. 

The  S-Curve  bornblet  is  comprised  of  an  axially  symmetric  body  or 
body-fin  configuration  which  is  designed  to  provide  an  unstable  restoring 
moment  at  small  angi e  of  attack  and  a  stable  pitching  moment  slope  at  a  large 
trim  angle  of  attack.  1'he  , S-Curve  name  is  derived  from  the  shape  of  the 
nonlinear  pitching  moment  curve.  / 

The  roll  -  through- zero  bornblet  is  comprised  of  a  lifting-body  configu¬ 
ration  equipped  with  a  roll  producing  device  which  will  cause  the  bornblet  to 
roll  in  a  direction  opposite  to  its  initial  direction  of  roll.  The  dispersion  of 
this  type  bornblet  is  inversely  proportional  to  the  magnitude  of  the  roll  torque 
used  to  reverse  the  roll. 

Doth  the  S-Curve  and  rid  1  -  through- ze  ro  bomblets  have  quasi  -balli  Stic 
flight  characteristics,  relatively  low  drag,  and  nose-first  impact. 

Preliminary  analytical  studies  of  the  flight  dynamics  and  dispersion 
characteristics  of  these  bomblets  were  accomplished  under  Air  Force  Con¬ 
tract  No.  F0863 5 -70 - C  - 00  J  2  and  the  results  presented  in  References  1  and 
2.  These  investigations  revealed  the  general  feasibility  of  both  the  S-Curve 
and  roi  1-th  rough  -  ze  ro  type  bomblets  and  provided  preliminary  estimates  of 
the  area  coverage. 

\ 

The  remaining  problems  re  (1)  to  determine  if  these  bomblets 
could  properly  function  under  representative  cluster  break-up  conditions, 

(2)  to  determine  the  performance  degradation  due  to  realistic  configurational 
and  aerodynamic  asymmetries,  (3)  to  determine  the  actual  impact  pattern 
distribution;  and  (4)  to  evolve  bornblet  shapes  and  candidate  configurations 
adaptable  to  a  wide  range  of  possible  warhead  and  weapon  Bystem  concepts. 


I 


The  present  effort  has  encompassed  both  a  wide  range  of  analytical 
and  experimental  programs  directed  toward  the  answers  to  these  problems. 
The  analytical  effort  has  entailed  the  development  of  a  comprehensive  Monte 
Carlo  trajectory  andampact  pattern  simulation  program,  based  on  modifica¬ 
tions  to  the  six-degrees-of-freedom  (6-DOF)  trajectory  program  (References 
3  and  4).  Concurrent  vKth  the  previous  and  present  analytical  efforts,  sever¬ 
al  series  of  wind  tunnel  tests  of  S-Curve  and  roll-through-zero  bomblet 
models  were  accomplished  by  the  Air  Force,  using  the  Arnold  Engineering 
Development  Center  (AEDC)  facilities.  These  data  have  been  integrated  into 
all  of  the  analytical  simulations  and  are  discussed  in  detail  in  the  body  of 
this  report. 

To  provide  additional  verification  of  the  S-Curve  bomblet  perform¬ 
ance,  a  flight  test  program,  using  gun-launched  models,  was  initiated  and 
completed  during  the  contractual  period.  The  test  plan  and  model  designs 
were  provided  by  the  contractor,  while  the  tests  and  data  acquisition  were 
accomplished  by  the  Air  Force  Armament  Laboratory.  A  summary  of  the 
test  results  and  analyses  of  these  data  are  provided  in  this  document. 

Because  the  aerodynamic  characteristics  of  bomblet  configurations 
undergoing  combined  rolling  and  coning  motions  at  large  trim  angles  of 
attack  are  still  not  well  understood,  the  development  of  a  special  dynamic 
wind  tunnel  test  system  was  undertaken.  The  test  system  was  designed  by 
the  contractor,  fabricated  by  the  Air  Force  Armament  Laboratory,  and 
installed  in  the  Armament  Laboratory  subsonic  wind  tunnel.  A  description 
of  the  test  system  and  its  operation  is  provided. 


SECTION  U. 


S-CURVE  BOMRLET  INVESTIGATIONS 

\  '  .  .  -r- . 

A.  CONFIGURATIONS 

The  present  effort  has  been  concerned  primarily  with  two  basic 
S-Curve  bomblet  configurations  as  depicted  in  Figure  1.  These  two  config¬ 
urations  were  selected  from  the  results  of  the  original  analytical  study 
(Reference  1)  as  having  the  best  performance  characteristics  of  those 
models  originally  evaluated  in  the  transonic  wind  tunnel  tests  (Reference  5). 

During  the  course  of  the  present  contractual  effort  some  additional 
configurational  modifications  to  the  basic  4-caliber  body  were  investigated 
in  the  wind  tunnel  facilities  of  AEDC,  including  two  different  blunt-nose 
shapes  and  fins  with  increased  span.  The  larger  fins  (0.14  caliber  exposed 
semi -span)  were  incorporated  as  a  means  of  improving  the  allowable  axial 
eg  range. 

Because  increased  span  fins  complicate  the  bomblet- dispenser  pack¬ 
aging,  design  studies  were  accompli  shed  to  determine  a  more  optimum 
shape.  Figure  2  illustrates  a  modified  afterbody  which  allows  the  exposed 
fin  semi -span  to  be  as  large*  as  -0.20  body-diameter  without  loss  of  packag¬ 
ing  efficiency.  Design  concepts  were  also  devised  for  rearward  extension 
of  the  fin  assembly  as  a  means  of  improving  the  packaging  arul  allowable 
axial  eg  range.  An  extensible  fin  design  is  depicted  in  Figure*  3.  The  con¬ 
figurations  shown  in  Figures  2  and  3  were  not  wind  tunnel  tested,  but  tests 
of  a  similar  extensible  tin  concept  for  an  S-Curve  version  of  the  BLU-87/B 
are  described  in  Reference  (>. 

Low  fineness  ratio  S-Curve  configurations  were  also  briefly 'exam - 
ined.  Figure  4  illustrates  two  low  fineness  ratio  shapes  with  S-Curve 
moment  characteristics  based  on  the  data  of  References  7  and  8.  The  lift, 
damping,  and  trim  stability  gradients  of  these  configurations  are  poor  com- 
'  pared  with  the  basic  4-caliber  designs.  ^ 

B.  AERODYNAMIC  CHARACTERISTICS  -  BASIC  S-CURVE 

CONFIGURATIONS 

The  original  aerodynamic  data  packages  for  the  basic  S-Curve  bomb- 
let  configurations  (Reference  1)  have  been  improved  and  extended  in  Mach 
number  as  a  result  of  further  wind  tunnel  tests  and  analyses.  This  work  is 
reviewed  in  the  following  paragraphs. 

Wind  Tunnel  Tests  The  static  and  dynamic  stability  characteristics 
of  the  two  basic  S-Curve  bomblet  configurations,  at  Mach  numbers  from  0.2 
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Figure*  I.  Basic  S-Curve  I’ki'rtblot  Configuration# 
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to  1.2,  were  previously  described  in  References  1,  5,  and  9.  In  this 
report,  only  the  results  of  more  recent  testing  will  be  described.  All  wind 
tunnel  tests  have  been  accomplished  at  the  Air  Force  Arnold  Engineering 
Development  Center  (AEDC)  test  facilities. 

During  the  period  of  the  present  contractual  effort,  additional  static 
and  dynamic  stability  data  for  the  basic  S-Curve  configurations  were  obtained 
at  Mach  numbers  1  .5,  2.0,  and  2.5.  Also,  magnus  data  and  fin  cant  effec¬ 
tiveness  were  obtained  at  subsonic  and  transonic  conditions  and  roll  damping 
data  were  measured  at  subsonic,  transonic,  and  supersonic  Mach  numbers. 
These  additional  series  of  tests  are  described  in  References  10  and  11. 

Subsequently,  two  new  blunt-nose  shapes  and  increased- span  fins 
were  evaluated  in  conjunction  with  both  the  cylindrical  and  boattail  afterbody 
models  (References  12  and  13).  Static  and  dynamic  stability  characteristics 
of  the  modified  configurations  were  determined  at  subsonic  and  transonic 
conditions. 

,  Static  Aerodynamic  Characteristics  The  supersonic  wind  tunnel 
data  for  the  two  basic  S-Curve  configurations  show  that  the  static  stability 
increases  at  supersonic  Mach  numbers.  The  S-Curve  moment  character¬ 
istics  are  retained  at  supersonic  Mach  numbers,  but  the  trim  angles  are 
less  than  those  at  subsonic  Mach  numbers.  The  trim  characteristics  of  the 
two  basic  bomblet  configurations  are  shown  in  Figure  5.  Note  that  the  trim 
data  for  the  boattail  bomblet  are  based  on  a  eg  position  0.228  calibers  aft 
of  that  for  the  cylindrical  afterbody  bomblet.  It  is  significant  that  the  trim 
angle  of  the  boattail  bomblet  varies  only  slightly  as  a  function  of  Mach 
number,  whereas  the  cylindrical  afterbody  configuration  experiences  large 
trim  changes  through  the  transonic  range. 

Dynamic  Stability  The  pitch  damping,  +  Cj^^. ,  increases  at 

supersonic  Mach  numbers,  for  both  basic  configurations.  There  is  only  a 
slight  effect  of  angle  of  attack  on  the  damping  at  supersonic  Mach  numbers. 

Magnus  Characteristics  The  magnus  tests  of  the  basic  cylindrical- 
afterbody  configuration  show  that  the  magnus  force  and  moment  are  consid¬ 
erably  less  than  previous  estimates  based  on  the  data  of  Reference  14. 

Magnus  force  and  moment  data  at  Mach  numbers  0.5  and  0.8  are  shown  in 
Figure  6.  The  data  were  obtained  for  a  test  Reynolds  number  of  550,  006 
based  on  body  diameter,  which  corresponds  to  a  flight  velocity  of  about 
500  feet/second  for  a  2-inch  diameter  bomblet  under  sea  level  conditions. 
Tests  at  various  Reynolds  numbers  show  that  the  magnus  data  are  sensitive 
to  Reynolds  number,  particularly  for  angles  of  attack  greater  than  about  16 
degrees. 

The  magnus  force  data  for  angles  of  attack  less  than  20  degrees  are 
positive  (negative  Cn^)  and  thus  are  in  agreement  with  body-alone  cross 
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Figure  6.  Magnus  Force  and  Moment  Characteristic s  of  Basic 
S-Curve  Bomblet  with  Cylindrical  Afterbody 
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flow  predictions  based  on  the  applicable  cross  flow  Reynolds  number.  For 
angles  of  attack  greater  than  eight  degrees  there  appears  to  be  some  magnus 
effect  due  to  fin-wake  interference,  since  the  magnus  force  is  different  for 
increasing  and  decreasing  spin  rates,  which  has  the  effect  of  reversing  the 
magnus  force  due  to  fin-wake  interference.  The  effect  of  nose  shape  on 
the  magnus  characteristics  is  very  pronounced  at  0.5  Mach  number,  and  at 
very  large  angles  of  attack  the  spherical  nose  bomblet  experiences  negative 
magnus  forces  at  Rj}^  =  550,000.  However,  the  magnus  force  becomes 
positive  if  the  Reynolds  number  is  decreased.  These  effects  are  believed 
to  be  due  to  the  fins,  since  the  body- alone  magnus  characteristics  should 
not  be  influenced  by  Reynolds  number  for  these  conditions. 

It  is  noteworthy  that  the  magnus  force  and  moment  are  quite  small  at 
the  nominal  trim  angles  of  attack  and,  hence,  the  tendency  for  coning  motion 
is  greatly  reduced  for  the  basic  S-Curve  configurations. 

The  reason  for  the  large  difference  between  the  present  results  and 
the  Australian  Weapons  Research  Establishment  data  on  a  similar  shape 
(Reference  14)  is  not  understood. 


The  magnus  characteristics  of  the  basic  S-Curve  configurations  at 
supersonic  Mach  numbers  were  estimated  from  data  correlations  on  similar 
shapes  and  from  consideration  of  both  body-fin  interference  and  cross-fk^•'• 
effects.  The  resulting  coefficients  were  found  to  be  highly  nonlinear  wi.' 
angle  of  attack. 


Spin  Characteristics  The  spin  damping  coefficient,  C,  ,  was 

P 

measured  for  the  basic  S-Curve  configuration  (Reference  11)  and  found  to 
be  less  than  the  original  estimates  presented  in  Reference  1,  although  the 
increase  in  damping  with  angle  of  attack  was  similar.  New  estimates  for 
the  spin  damping  at  supersonic  velocities  Mere  required,  since  spin  damping 
tests  were  not  accomplished  for  Mach  numbers  greater  than  0.9.  The 
supersonic  damping  was  estimated  from  the  subsonic  data,  using  propor¬ 
tionality  factors  based  on  the  fin  normal  force  derivative. 


The  fin  effectiveness  coefficient,  ,  was  measured  late  in  the 

program  for  the  model  with  large  span  fins  (Reference  13).  These  data  are 
shown  in  Figure  7.  These  coefficients,  when  corrected  for  the  smaller 
siae  basic  fins,  are  still  somewhat  larger  than  the  estimated  curve,  which 
was  used  for  most  of  the  motion  simulations. 


Effect  of  Aerodynamic  H<.  !’  Angle  During  the  present  effort,  the 
effect  of  aerodynamic  roll  angle  on  configuration  Fg,  was  inves¬ 

tigated  in  the  wind  tunnel  at  supersonic  Mach  numbers  (Reference  10).  The 
first  harmonic  coefficients,  Cgjrj,  and  C  ^  ^  ,  were  evaluated  at 

s  1 5  degrees.  For  angles  of  attack  less  than  10  degrees  there  is  no 
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significant  effect  due  to  aerodynamic  roll  angle,  but  at  larger  angles  of 
attack  these  coefficients  increase  rapidly  and  become  very  nonlinear  with 
angle  of  attack.  The  induced  roil  and  side  moment  coefficients,  CSM  and 
Cfc$j  »  are  Pitted  in  Figure  8  and  compared  *»ith  the  subsonic  results. 

Effect  of  Asymmetrical  No.se-Roughneas  It  has  been  reported  that 
asymmetrical  nose  roughness  can  have  a  significant  effect  on  the  side  force 
side  moment,  and  roll  moment  (Reference  15),  If  nose  roughness  is  confined 
to  a  single  meridional  orientation,  the  above  coefficients  can  experience  an 
abrupt  change  in  magnitude  and  sign  when  the  body-fixed  roughness  plane  and 
the  angle  of  attack  plane  coincide.  A  typical  variation  of  the  side  moment 
coefficient  variation  as  a  function  of  the  aerodynamic  roll  angle  is  depicted 
tn  Figure  9.  r 


To  model  such  affects  in  the  6-DOF  computer  program,  a  special 
aerodynamic  step  function  was  provided,  which  together  with  first  and  second 
order  roll  dependent  harmonics,  permitted  a  close  approximation  of  the  test 
data.  Figure  9  also  shows  a  typical  fit  obtained  with  the  step  function  and 
harmonics  Provision  was  also  made  for  the  step  function  to  have  a  differ¬ 
ent  phase  for  the  aide  moment  and  the  roll  moment. 


The  asymmetrical  nose- rough 
of  the  6-DOF  motion  simulations,  but 
disiMict  aerodynamic  asymmetry. 


ness  effect  was  not  incorporated  for  all 
was  investigated  as  a  separate  and 


°f  Lqsisaassf.^./’in*  The  effect  of  increasing  the  c-xuosed 
sem:  -  span  ot  ih«  basic  S  Curve  bombiet  fins  from  0.078  caliber  to  0  H 

cahbur  «,  inv,«ls«,d  a.  »  mo**.  ...  ,mprov,ns  th.  bomblw  allowable  «, 
rang,,  lb.  oili.ct  ol  th.  larji.r  I'lns  on  'hi-  normal  Ion  .,  center  o.  pressure 
is  depicted  in  Figure  10.  The  increased  ho  span  allows  a  0.1  caliber  aft 
cf  for  the  cylindrical  afterbody  configuration  ami  a  0.15  caliber  aft  ce 
shift  for  the  tmattaU  afterbody  configuration  at  subsonic  conditions.  The 
larger  fms  allow  the  eg  to  be  .at  or  near  the  body  midpoint.  Data  are  also 
shown  for  a  configuration  with  im  leased  nose  curvature,  Although  the 
center- ef-pressu  re  is  moved  rearward,  the  S  Curve  moment  character- 

t*tic  is  lost,  since  the  center  of  pressure  moves  forward  with  increasing 

angle  of  attack.  8 


The  increase  in  fin  span  was  found  to  have  negligible  effect  on  the 
induced  aerodynamic  forces  and  moment  coefficients,  Cc,- .  (!al  ^nd 

rh  i!rjAI18  CS  att*ck  a  degrees.  For  angles  of  attack  larger 

than  20  degrees  there  is  some  difiVtence  in  the  induced  aerodynamic  coeffi¬ 
cients  for  the  two  different  fin  sixes,  particularly  the  roll  moment.  C* 

nr 

The  magnus  force  coefficients  for  configuration  Rs%2ASl  .  with 
the  increased  span  fins,  are  about  equal  or  less  than  those  for  the  smaller 
fins  at  angles  of  attack  less  than  eight  degrees.  The  Magnus  data  for 


Figure  %  fcffeet  of  Nosr-KeughneftiJ  the  tmiure<£  ArnHynimte 
Jsttie  Moment 
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configuration  BgNg2AgjFg5  are  described  in  Reference  13.  At  angles  of 
attack  greater  than  twenty  degrees  the  magnus  force  becomes  increasingly 
positive  (C jsj  <  0)  for  the  large  fins.  The  magnus  moment,  in  contrast,  is 

significantly  increased  for  the  larger  fins  at  all  angles  of  attack  less  than 
about  twenty  degrees,  and  Cjy{p  is  positive  throughout  this  angle  of  attack 

range  for  Mach  numbers  of  0.9  and  less.  The  magnus  moment  coefficients 
for  the  configuration  with  large  fins  are  highly  nonlinear  with  angle  of  attack, 
and  CMp  becomes  negative  for  angles  of  attack  greater  than  twenty  degrees. 

Interestingly,  positive  fin  cant  reduces  the  magnus  moment,  and  for  3  de¬ 
grees  cant  the  magnus  moment  is  near  zero  for  subsonic  Mach  numbers  and 
angles  of  attack  less  than  about  10  degrees. 

Aerodynamic  Coefficients  -  6-DOF  Simulations  A  tabulation  of  the 
aerodynamic  coefficients  used  for  the  6-DOF  motion  simulations  are  pre¬ 
sented  in  Tables  I  and  II  for  the  basic  cylindrical  afterbody  and  boattail 
afterbody  configurations,  respectively.  The  aerodynamic  coefficients  are 
tabulated  both  as  a  function  of  Mach  number  and  angle  of  attack.  Coefficients 
are  shown  for  only  five  Mach  numbers  because  of  the  limitation  of  the  6-DOF 
trajectory  program.  A  detailed  description  of  the  coefficient  definitions  can 
be  found  in  References  3  and  16.  The  coefficients  are  divided  into  two 
groups;  the  first  group  of  coefficients  represents  the  aerobalii stic  coeffi¬ 
cients  which  are  independent  of  the  y-z  lx>dy~axis  roll  orientation.  The 
second  group  of  coefficients  is  applicable  only  to  a  specific  set  of  y-z  body- 
fixed  axes,  which  are  aligned  with  a  plane  of  rotational  symmetry.  This 
second  group  of  coefficients  defines  the  induced  aerodynamic  forces  and 
moments,  which  are  harmonic  functions  of  the  aerodynamic  roll  angle,  as 
well  as  the  coefficients  which  describe  the  configurational  and  aerodynamic 
asymmetries. 


C.  PREDICTED  FRIGHT  DYNAMICS  OF  BASIC  CONFIGURATIONS 

Effect  of  Roll  Dependent  Aerodynamic  Coefficients  The  addition  of 
roll  dependent  aerodynamic  coefficients,  Cjjj,-,  and  Cy.^,  to  the  pre¬ 

vious  6-  IX) F  motion  simulations  (Keteronce  1)  produced  some  interesting 
results.  *  For  small  configurational  asymmetries  and  small  cant  angles, 
it  was  fo  nd  that  roll  lock-in  was  occurring  for  nearly  all  of  the  computed 
motion  histories.  It  was  also  noted  that  the  lock-in  was  not  a  direct  result 
of  the  roll  torques.  It  was  discovered  that  the  lock-in  is  critically  affected 
by  the  roll  dependent  side  moment.  A  re-examination  of  the  equation  for 
zero  coning  motion  (Appendix  I)  has  shown  that  for  zero  roll  rate  the  roll 
dependent  side  moment  must  have  positive  gradient  with  respect  to  the 


1 


The  induced  aerodynamic  coefficients  were  not  included  in  the  original 
S-Curve  bomblet  simulations,  Reference  1. 
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TABLE  II.  AERODYNAMIC  COEFFICIENT  SUMMARY  FOR  BASIC 
S-CURVE  CONFIGURATION  BsNslAS2FS3 
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aerodynamic  roll  angle  for  the  zero  coning  solutions  to  be  stable,  i.  e.  , 


— — —  >  0  :  stable  zero  coning 

Figure  li  illustrates  the  side  moment  coefficient  variation  with  4>  for  con¬ 
figuration  bsnS2as1fs2  and  the  values  of  $  at  lock-in  for  a  series  of  ten 
trajectories  which  included  both  roll  dependent  moments  and  eg  lateral 
offset.  In  all  cases  lock-in  occurs  in  regious  where  the  inequality  is  satis¬ 
fied. 

The  lock-in  condition  implies  that  the  roll  rate  and  coning  rates  are 
identical,  so  that  zero  coning  is  possible  only  for  the  singular  case  of  zero 
roll  rate.  Although  the  side  moment  controls  the  region  of  lock-in,  the 
specific  lock-in  conditions  are  determined  by  the  roll  torques.  The  effect 
of  the  roll  torques  will  be  discussed  next. 

Effect  of  CG  Lateral  Offset  Since  for  small  cant  angles  lock-in 
can  be  expected  from  the  side  moment,  as  described  above,  it  is  important 
to  see  how  die  lock-in  may  influence  the  roll  dynamics,  since  the  roll  rate 
and  coning  rate  will  be  identical.  The  roll  dynamics  can  be  greatly  influ¬ 
enced  by  eg  lateral  offset,  as  well  as  by  the  fin  cant  and  fin-induced 
aerodynamic  moments.  From  the  diagram  of  Figure  12,  it  is  seen  that 
the  roll  moment  due  to  eg  offset  is^ 

q  '  CN  sin  f  (1) 

where 

f,  -  it/2  -  $  *  c. 

Thus,  the  roll  moment  due  to  eg  offset  depends  both  upon  the  aerodynamic 
roll  angle,  4*  ,  and  the  orientation  of  the  fin  symmetry  planes  with  respect 
to  the  orientation  of  the  eg  offset,  which  is  given  by  the  am  ’<  *  ,  ppr  a 

particular  bomblet  <;  is  a  constant,  while  the  aerodynamic  roil  angle,  $  , 
is  influenced  strongly  by  the  fact  that  the  motion  tends  toward  a  condition 
where  d  Cyq^/d  *  has  a  positive  gradient.  This  implies  that  the  angle  ?,  , 
which  controls  the  roll  moment  due  to  eg  lateral  offset,  wilt  have  discon¬ 
tinuous  ranges  of  values,  and  consequently,  the  induced  roll  moment  itself 
will  change  incrementally,  depending  upon  the  lock-in  angle.  The  induced 
roll  moment  due  to  eg  offset  reaches  a  maximum  when  f  has  values  of 
tt/2,  W2,  etc. 

> 

Equation  (If  assumes  that  all  other  transverse  aerodynamic  forces  are 
zero. 


ornent  on  Roll  Lock-In 


"  '"J.'  •  v'" gsis»TO**F! ; 

»»*iniii'f  r 


Figure*  12.  Definition  of  Angles  and  Axes  Used  in  Describing 
Bomblet  Roll  Dynamics 


The  roll  rate  (and  also  the  coning  rate)  will  obviously  depend  upon 
both  the  roll  driving  torques  and  roll  damping,  although  the  direction  of  roll 
will  in  many  instances  be  determined  by  the  roll  torque  due  to  eg  offset, 
because  it  will  often  be  larger  than  the  fin  cant  and  fin  induce*1  aerodynamic 
roll  moments. 


Detailed  examination  of  6 - IX) F  motion  histories  for  a  series  of  runs 
where  the  standard  deviation  of  &y  was  0,008  inch,  showed  that  for  all  the 
runs  where  lock-in  occurred,  the  tetminal  spin  rate  was  of  the  same  sign 
and  approximately  proportional  to  the  induced  roll  moment  resulting  from 
the  eg  lateral  offset. 

The  effect  of  eg  lateral  offset  on  the  dispersion  is  shown  in  Figure 
1 1.  Th«*  data  were  obtained  from  6-DOF  trajectory  calculations  with  and 
without  eg  offset ,  and  show  the  increase  or  decrease  in  radial  dispersion 


C.G.  OFFSET  (Inchtt) 


DECREASED  - -  - *•  INCREASED 

DISPERSION  DISPERSION 

CHANGE  IN  RADIAL  DISPERSION  (F«tt) 

Figure  13.  Kffeet  of  CG  Lateral  Off M>t  on  Radial  Dispersion 


as  a  function  of  eg  lateral  offset.  As  can  he  seen,  the  dispersion  data 
tend  to  he  erratic  because  of  the  incremental  changes  in  the  roll  lock-in 
angle  and  roll  rate.  It  is  also  noteworthy  that  a  eg  offset  of  only  0.002 
inch  has  a  large  effect  on  the  dispersion. 

Effect  of  Fin  Cant  For  large  eg  offsets,  the  preceding  indicates 
that  large  roll  anti  coning  rates  can  be  anticipated,  as  long  as  lock-in  occurs. 
An  obvious  approach  to  the  prevention  of  lock-in  is  the  use  of  large  intention¬ 
al  fin  cant,  such  that  the  effects  of  eg  offset  and  the  side  moment  are  over¬ 
come. 

A  preliminary  evaluation  was  made  using  a  cant  angle  of  0.5  degree 
an«l  a  standard  deviation  of  eg  lateral  offset  of  0.008  inch.  For  a  trim 
angle  of  attack  of  10  degrees  and  a  eg  offset  of  0.008  inch,  the  maximum 
roll  torque  due  to  eg  offset  is  approximately  the  same  as  that  due  to  cant. 


2.3 


Out  of  10  flight  simulations  only  one  lock-in  condition  was  noted.  Figure  14 
compares  the  roll  rates  and  angle  of  attack  plane  histories  for  two  bomblet 
flights  differing  only  in  the  cant  angle  and  eg  offset.  The  angle  of  attack 
plane  rotation  (coning)  is  significantly  reduced  with  the  intentional  fin  cant, 
even  though  the  roll  rate  is  greatly  increased. 

The  effectiveness  of  intentional  fin  cant  has  also  been  proven  from 
large  Monte  Carlo  samples  where,  in  addition  to  eg  lateral  offset,  lateral 
misalignment  and  nose -roughness  asymmetry  were  also  present. 

It  should  be  mentioned  that  the  effectiveness  of  intentional  fin  cant  is 
due  in  part  to  the  relatively  small  spin  dependent  magnus  moment.  This 
allows  the  roll  rate  to  be  Increased  to  moderate  values  without  adverse 
magnus  effect  on  the  coning  motion. 

The  maximum  cant  angle  is  established  either  by  magnus  considera¬ 
tions  or  by  the  possibility  of  the  configuration  becoming  gyroscopically 
stable.  The  cant  angle  required  for  gyroscopic  stability  decreases  very 
rapidly  as  the  bomblet  fineness  ratio  is  decreased,  due  both  to  the  decrease 
in  the  ratio  of  the  transverse  to  axial  moment  of  inertia  and  the  reduction  in 
unstable  pitching  moment  derivative  at  zero  angle  of  attack. 

D.  DISPERSION  PREDICTION  -  GENERAL  CONSIDERATIONS 

It  is  difficult  to  assess  the  impact  dispersion  characteristics  of  the 
S- Curve  type  bomblet  without  consideration  of  the  exact  bomblet  motion. 
This  is  due  to  the  fact  that  the  aerodynamic  lift  force  (which  provides  the 
dispersion)  is  constantly  changing  in  magnitude  and  direction  as  a  result  of 
the  bomblet  dynamics.  The  dynamics  of  the  S- Curve  type  bomblet,  in  turn, 
are  highly  sensitive  to  the  initial  flight  conditions  and  bomblet  configura¬ 
tional  asymmetries. 

For  these  reasons  the  trajectory  calculations  have  been  formulated 
so  as  to  provide  the  closest  possible  duplication  of  actual  bomblet  flight. 

This  has  necessitated  that  the  initial  flight  conditions  at  dispi  user  opening 
and  during  cluster  break-up  be  realistically  simulated,  and  that  bomblet 
configurational  asymmetries  be  introduced  which  are  representative  of 
actual  production  ordnance. 

In  contrast,  a  simplified  treatment  of  the  bomblet  dispersion  based 
on  constant  trim  angle  of  attack  and  no  angular  motion,  would  not  only  over¬ 
predict  the  dispersion  but  would  not  provide  an  impact  pattern  distribution. 
Subsequently,  poor  agreement  between  the  simulations  and  flight  results 
would  be  expected. 
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CG  OFFSET  FIN  CANT 

-  NONE  0.05  DEGREES  MISALIGNMENT 

-  0.011  INCHES  0.50  DEGREES  (INTENTIONAL) 

V  (Dtgrtct) 


p  ( Rodiont/Stcond) 


Figure  14.  Kffect  of  Intentional  Fin  Can!  on  the  Motion  of  an 
S-Curve  Hombict  with  Lateral  CG  Offset 


Monte  Carlo  Method  To  the  maximum  practicable  extent,  the  per¬ 
formance  of  the  S-Curve  bomblet  has  been  investigated  using  die  Monte 
Carlo  method,  taking  into  account  both  the  randomness  of  the  initial  flight 
conditions  and  the  configurational  asymmetries.  The  Monte  Carlo  technique 
and  its  adaptation  to  the  simulation  of  bomblet  cluster  break-up  and  configu¬ 
rational  asvmmetries  is  described  in  Section  IV.  Except  for  center-of- 
gravity  lateral  offset,  parametric  analyses  (with  one  or  more  fixed  values 
of  each  parameter)  have  been  avoided,  since  results  obtained  in  this  way 
can  be  misleading  unless  a  sufficient  range  of  all  variables  is  investigated. 
The  limitation  of  the  parametric  type  of  analysis  is  simply  a  result  of  the 
fact  that  nearly  all  of  the  parameter  responses  are  nonlinear.  This  pre¬ 
cludes  any  generalization  of  the  results  or  the  use  of  linear  superposition 
as  a  means  of  determining  the  combined  effect  of  a  large  number  of  param¬ 
eters.  Although  parametric  analysis  is  widely  used  in  the  determination  of 
the  dispersion  of  conventional  bombs,  projectiles,  and  missiles,  it  is  not 
appropriate  for  self-dispersing  ordnance. 

E.  MONTE  CARLO  IMPACT  PATTERNS 

Except  where  noted,  dispersion  data  have  been  computed  for  bomb- 
lets  of  standardized  size,  weight,  and  inertia.  The  standardized  physical 
characteristics  data  are  summarized  in  Table  III.  Also,  emphasis  has 
beenpi  vcedon  three  dispenser  event  conditions,  corresponding  to  dispenser 
or  warhead  function  at  high  subsonic,  transonic,  and  supersonic  Mach 
numbers,  -espectively.  The  high  subsonic  event  condition  is  defined  as 
the  basic  event  condition,  and  is  identical  to  that  considered  for  the  original 
S-Gurve  dispersion  studies  (Reference  1).  The  transonic  event  conditions 
are  representative  of  those  for  a  prototype  guided  dispenser,  while  the 
supersonic  event  conditions  are  based  on  the  flight  characteristics  of  a 
surface-to-surface  tactieal  missile  warhead. 


TABLE  III.  PHYSICAL  CHARACTERISTICS  DATA  FOR 
STANDARD  SI ;*E  S- CURVE  BOMBLET 

Length 

8  inches 

Diameter 

2  inches 

Aerodynamic  Reference  Area 

0.0218 

Aerodynamic  Reference  Length 

0.  lb?  feet 

Weight 

I.^I  pounds. 

Axial  Moment  ot  Inertia,  !x 

Transverse  Moment  of  Inertia.  I 

l.A2d  x  10‘4  slug -ft2  1 

1 . Of.  x  !0'*^  slug* st" 

Transverse  Moment  of  Inertia.  I 


Except  where  noted,  the  bomblet  initial  motion  perturbations  at 
cluster  break-out  and  the  bomblet  configurational  asymmetries  correspond 
to  the  values  described  in  Section  IV  of  this  report. 

Impact  Dispersion  Pattern  -  Basic  Event  Condition  Figure  1  5 
shows  the  computed  Monte  Carlo  impact  pattern  data  for  the  basic  high- 
subsonic  event  condition,  with  2000  feet  height-of-burst  (HOB).  The  pattern 
is  computed  for  the  cylindrical  afterbody  bomblet  configuration 
®SNS2^S1  Fs2*  with  the  center-of- gravity  positioned  for  a  nominal  subsonic 
trim  angle  of  attack  of  10  degrees.  The  impact  pattern  is  a  composite  of 
two  separate  simulations.  The  open  symbols  represent  a  Monte  Carlo  sim¬ 
ulation  of  100  runs  in  which  the  eg  lateral  offset  was  neglected,  but  all 
other  configuration  asymmetries  were  included.  The  solid  symbols  repre- 
sent  a  simulation  of  70  runs  with  a  standard  deviation  of  eg  lateral  offset 
of  0.0016  inch.  In  the  latter  series  the  fin  cant  was  identically  zero.  The 
standard  deviation  values  of  cross -range  dispersion  are  within  two  percent 
for  the  two  groups,  while  the  standard  deviation  values  of  range  dispersion 
are  within  ten  percent  for  the  two  groups.  Thus,  small  values  of  eg  lateral 
offset  do  not  significantly  affect  the  impact  pattern,  even  though  the  effect  on 
individual  trajectories  may  be  large. 

The  overall  pattern  length  is  seen  to  be  about  1000  feet,  while  the 
overall  pattern  width  io  about  800  feet.  For  the  combined  patterns,  the 
standard  deviations  of  range  and  cross- range  dispersion  with  respect  to  the 
mean  center  of  impact  are  231  ami  183  feet,  respectively.  This  shows  that 
the  impact  points  tend  to  have  greater  concentration  toward  the  center  of  the 
pattern.  Figure  16  shows  the  cumulative  probability  distributions  for  the 
first  group.  The  distribution  of  impact  points  is  piecewise  normal,  ^  em¬ 
phasizing  the  fact  that  the  impact  dispersion  characteristics  are  the  result 
of  nonlinear  interactions  between  the  dispersion  parameters. 

Comparison  with  /, e r o  Coning  Solutions  It  is  of  interest  to  compare 
the  Monte  Carlo  results  with  a  pattern  prediction  based  on  zero  coning.  For 
the  basic  event  condition,  a  scries  of  zero  coning  trajectories  were  com¬ 
puted  with  the  initial  orientation  of  the  angle  of  attack  plane  as  a  parameter. 
The  initial  angle  of  attack  was  in  all  cases  assumed  to  be  *>  degrees  (nose 
outward  from  the  line  of  flight)  and  cluster  break-up  was  assumed  to  occur 
instantaneously  at  an  altitude  200  feet  below  the  dispenser  HOB.  This  point 
represents  approximately  the  lowest  cluster  break-out  point  in  the  Monte 
Carlo  simulations.  The  locus  of  impact  points  computed  in  this  manner  is 
shown  superimposed  on  the  Monte  Caito  pattern  of  Figure  16. 


5  The  normal  distribution  appears  as  a  straight  line  on  the  probability 
graph. 
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PERCENT  CUMULATIVE  PROBABILITY 


It  is  apparent  that  the  simplified  theory  overpredicts  the  Monte  Carlo 
dispersion  pattern;  however,  the  maximum  dispersion  values  computed  in 
the  Monte  Carlo  simulations  closely  approach  the  zero  coning  boundary. 

Transonic  and  Supersonic  Delivery  Monte  Carlo  dispersion  data 
for  the  transonic  dispenser  event  condition  were  computed  for  the  basic 
boattaii-aiLerbody  bomblet,  because  this  configuration  has  much  better  trim 
characteristics  through  the  tfahsonic  range.  Preliminary  pattern  size  pre¬ 
dictions  based  on  zero  coning  showed  that  the  boattail  bomblet  would  have 
about  twice  the  dispersion  of  the  cylindrical  afterbody  bomblet.  The  cluster- 
breakout  perturbations  and  configurational  asymmetries  correspond  to  the 
data  of  Section  IV,  except  that  eg  lateral  offset  is  neglected. 

Figure  17  shows  the  computed  impact  pattern  for  the  transonic  dis¬ 
penser  event  conditions  with  HOB  of  2000  feet.  The  standard  deviations  of 
cross-range  and  range  dispersion  with  respect  to  the  mean  center  of  impact 
are  within  ten  percent  of  those  for  the  subsonic  event  condition. 

For  the  supersonic  event  (Vjr  =  2000  fps)  the  boattail  configuration 
was  again  selected.  The  impact  pattern  data  are  shown  in  Figure  18.  The 
results  show  a  very  large  increase  in  the  area  coverage,  with  a  few  bomb- 
lets  impacting  several  thousand  feet  forward  of  the  mean  center  of  impact. 
Thus,  it  is  evident  that  the  large  dynamic  pressure  can  be  exploited  for 
increased  dispersion,  notwithstanding  the  proportionately  larger  terminal 
spin  rates  and  tendencies  for  coning  motion. 

liable  IV  summarizes  *-he  statistical  dispersion  data  for  the  subsonic, 
transonic,  and  supersonic  event  conditions  which  have  been  evaluated.  It 
will  be  noted  that  the  standard  deviation  values  of  the  range  and  crOoS- range 
dispersion  for  the  supersonic  event  condition  are  greater  than  the  subsonic 
and  transonic  values  by  factors  of  about  four  and  two,  respectively.  Ob¬ 
viously,  this  larger  dispersion  can  also  be  exploited  by  a  decrease  in  the 
HOB. 

Effect  of  Static  Unbalance  The  static  unbalance  (lateral  center-of- 
gravity  offset)  of  a  bomblet  configuration  can  be  expected  to  vary  widely 
depending  upon  the  method  of  manufacture  and  assembly.  Using  the  sub¬ 
sonic  event  condition,  separate  Monte  Carlo  impact  pattern  simulations 
were  accomplished  for  the  following  values  of  static  unbalance:  0,  0.0016, 
0.0080,  and  0.0400  inch.  The  initial  motion  perturbations  and  lateral  mis¬ 
alignment  remained  identical  to  those  used  for  the  previously  described 
simulations.  A  minimum  of  70  trajectories  was  computed  for  each  simu¬ 
lation.  For  these  simulations  the  fin  cant  was  assumed  to  be  identically 
zero.  The  results  of  the  simulations  are  summarized  in  Figure  10,  which 
is  a  plot  of  both  the  maximum  dispersion  and  the  standard  deviation  of  dis¬ 
persion  versus  the  eg  lateral  offset  parameter.  The  results  show  that  eg 
offsets  greater  than  0.001  inch  have  a  measurable  effect  on  the  dispersion 
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TABLE  IV.  COMPARISON  OF  IMPACT  DISPERSION  DATA  FOR 
SUBSONIC,  TRANSONIC,  AND  SUPERSONIC  EVENT  CONDITIONS 


Event  Conditions 

. . .  ~ 

Mean 

Range 

X 

( Feet) 

Mean 

Cross -Range 

Y 

( Feet) 

Standard 

Deviation 

Range 

aX 

{ Feet) 

Standard 

Deviation 

Cross-Range 

°Y 
( Feet) 

Subsonic  Event 

VE  =  900  fps 

YE  =  45" 

1853 

.17 

221 

181 

Transonic  Event 

VE  =  114  5  fps 

YE  ■  65.7" 

87  5 

1  1 

202 

183 

Supersonic  Event 

VE  2000  fps 

.  Y  E  50" 

1978 

-M9 

959 

47  0 

Note:  l)  HOB  ,’.000  feet. 
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and  that  for  eg  offsets  greater  than  0.01  inch  a  very  large  decrease  in 
both  the  maximum  and  standard  deviation  of  dispersion  can  be  expected. 

Effect  of  Fin  Cant  Some  of  the  preceding  simulations  were  repeated 
using  intentional  values  of  fin  cant  of  0.1,  0.5,  and  2.5  degrees.  The  0.5 
degree  cant  angle  was  found  to  significantly  increase  the  dispersion  for  eg 
lateral  offsets  of  0.0016  and  0.0080  inch,  and  the  results  are  superimposed 
on  the  zero  cant  data  in  Figure  19.  The  0.1  degree  cant  angle  was  found  to 
be  insufficient  to  prevent  roll  lork-in  even  for  a  eg  offset  of  0.0016  inch, 
and  out  of  35  trajectories,  21  developed  a  roll  lock-in  type  motion  with 
constant.  The  2.5  degree  cant  angle,  which  was  investigated  in  conjunction 
with  the  0.0400  inch  eg  offset,  did  not  improve  the  dispersion.  This  was 
due  both  to  the  increased  magnus  moment  and  the  tendency  toward  spin 
stabilization. 

It  should  be  emphasized  that  a  finite  trim  angle  of  attack  due  to  the 
S-shaped  moment  curve  will  exist  only  if  the  configuration  is  gyroscopically 
unstable.  The  cant  angle  required  for  gyroscopic  stability  at  zero  angle  of 
attack  is  only  3.0  degrees  at  Mach  number  0.5.  Because  of  the  nonlinearity 
of  the  moment  curves,  the  non- rolling  trim  is  reduced  by  cant  angles  less 
than  3.  0  degrees. 

Eftect  ot  Initial  Pitch  Disturbance  For  tilt*  basic  Monte  Carlo  sim¬ 
ulations  an  initial  pitch  rate  at  cluster  break-out  was  introduced  such  that  the 
magnitude  of  the  first  peak  in  the  angle  of  attack  oscillation  had  a  2o  value 
ot  about  45  degrees.  For  a  900  fps  event  velocity,  the  corresponding 
standard  deviation  for  the  pitch  rate  was  estimated  to  be  about  3C>  radians/ 
second.  Statistical  analysis  of  30  Monte  Carlo  simulations  computed  with  a 
standard  d-viatmn  of  initial  pitch  rate  equal  to  35  radians / second  showed 
that  the  maximum  angle  of  attack  was  52  degrees  and  that  the  mean  value  of 
the  first  oscillation  peak  was  30  degrees.  A  second  Monte  Carlo  simulation 
was  accomplished  with  twice  the  pitch  rate  at  cluster  break-out.  Again, 
from  analysis  of  30  trajectories,  the  maximum  angle  of  attack  was  found  to 
be  9?  degrees.  hiie  the  on  value  of  the  first  oscillation  peak  was  in¬ 
creased  to  46  de,  rees  with  a  standard  deviation  of  21  degrees. 

The  effect  of  the  increased  pitch  rate  on  the  impact  dispersion  was 
surprisingly  small;  the  standard  deviation  of  cross-range  dispersion  was 
d<*ci  *a  <’d  About  20  percent,  while  the  standard  deviation  of  range  disper¬ 
sion  u.  i  i dually  increased  about  5  percent. 

ffect  of  Increased  Trim  Au.,lc  of  Attack  The  impact  patterns 
shown  Figures  15,  17.  and  18  arc  based  on  bomblet  configurations  which 
have  a  nominal  trim  angle  of  attack  of  approximately  10  degrees  at  low  sub¬ 
sonic  velocities.  This  trim  angle  was  selected  conservatively,  so  as  to 
avoid  possible  anomalous  bomblet  motions  due  to  the  highly  nonlinear  induced 
aerodynamic  moments  at  large  angles  of  attack.  To  determine  the 
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significance  of  these  nonlinearities  at  large  angles  of  attack,  the  nominal 
trim  angle  of  attack  of  the  cylindrical -afterbody  bomblet  was  increased  to 
16  degrees,  and  a  new  Monte  Carlo  simulation  accomplished.  The  impact 
patterns  are  compared  in  Figure  20.  For  the  large  *-rim  angle,  the  standard 
deviation  values  of  range  and  cross-range  dispersion  are  increased  by  fac¬ 
tors  of  1.67  and  2.20,  respectively.  These  compare  to  a  factor  of  1,72  for 
the  bomb.el  normal  force  increase  at  the  larger  trim  angle. 

Effect  of  Bomblet  Mass  Density  While  the  S-Curve  bomblet  disper¬ 
sion  is  essentially  proportional  to  the  lift/weigh::  ratio,  the  effect  of  the 
bomblet  mass  and  inertia  is  also  felt  through  the  pitch  and  roll  dynamics. 

To  determine  more  exactly  the  relationship  between  the  dispersion  and  the 
bomblet  mass  density,  Monte  Carlo  simulations  were  accomplished  for  two 
bomblet  configurations  differing  only  in  mass  density.  The  heavy  bomblet 
was  assumed  to  have  twice  the  mass  density  of  the  standard  bomblet.  For 
these  simulations  the  boattail  bomblet  configuration  was  assumed,  with  the 
nominal  trim  angle  of  attack  adjusted  to  sixteen  degrees.  The  basic  sub¬ 
sonic  event  condition  was  selected  for  the  analysis. 

The  computed  impact  patterns  (based  on  30  trajectories  each)  were 
statistically  analyzed  and  it  was  found  that  the  standard  deviations  of  range 
and  cross-range  dispersion  for  the  heavy  bomblet  were  47  percent  of  those 
for  the  reference  bomblet;  in  other  words,  the  dispersion  was  approximately 
proportional  to  the  reciprocal  of  the  mass  ratio. 

The  results  show  that  the  effects  of  mass  on  the  flight  time  and  flight 
velocity  tend  to  be  compensatory  insofar  as  the  dispersion  is  concerned. 
Also,  the  dispersion  is  apparently  not  very  sensitive  to  changes  in  the  roll 
and  pitch  dynamics  resulting  rrom  the  increased  moments  of  inertia. 

Effect  of  Nose- Koughnesa  Asymmetry  The  effect  which  nose- 
roughness  asymmetry  can  have  cm  the  aerodynamic  characteristics  of  an 
S-Curve  type  bomblet  is  d  scribed  in  Section  Ii-B.  To  determine  the  effect 
of  nose- roughness  type  asymmetry  on  the  dispersion,  Monte  Carlo  simula¬ 
tions  were  accomplished  for  the  cylindrical -afterbody  bomblet  configuration 
using  the  basic  subsonic  event  conditions.  The  standard  deviation  values  for 
cant  and  lateral  misalignment  were  identical  to  those  derived  in  .Section  IV, 
while  the  standard  deviation  of  eg  lateral  offset  was  assumed  to  be  0.0016 
inch. 

For  the  first  Monte  Carlo  simulation  the  mean  fin  cant  was  set  to 
zero.  The  impact  data,  in  the  case  of  zero  mean  cant,  showed  a  40  percent 
reduction  in  the  standard  deviation  values  of  range  and  cross- range,  when 
comparison  was  made  with  the  impact  pattern  data  computed  without  nose- 
royghness  asymmetry.  Detailed  examination  of  the  bomblet  motion  histories 
showed  that  periods  of  roll  lock-tt*  and  erratic  pitching  motion  were  encount¬ 
ered. 
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Figure  20.  Effect  of  S-Curvo  Homblet  Trim  Angle  of  Attack  on  Monte  Carlo 
Impact  Pattern  for  Standard  Subsonic  Event  Condition 
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The  Monte  Carlo  simulation  was  repeated  with  a  mean  fin  cant  angle 
of  0.5  degree.  With  intentional  fin  cant  the  reduction  in  the  standard  devia¬ 
tion  of  cross-range  dispersion  was  only  3  percent  while  the  reduction  in 
range  dispersion  was  but  15  percent. 

Thus,  nose-roughness  asymmetry  should  not  be  of  great  concern  if 
intentional  fin  cant  is  utilized. 

F.  SPECIAL  STUDIES  -  HIGH  DRAG  AND  DUAL  MODE  BOMBLETS 

For  certain  bomblet  warheads,  a  near-vertical  flight  path  at  impact 
is  desirable  as  a  means  of  achieving  optimum  weapon  effectiveness. 

High-Drag  Bomblet  A  brief  investigation  was  made  to  determine 
the  effect  of  drag  on  the  flight  characteristics  of  a  typical  S-Curve  bomblet. 
The  high-drag  configuration  was  assumed  to  have  a  two-caliber  attached- 
inflatable-decelerator  (AID),  which  would  be  configured  such  that  the  bomb- 
let  would  retain  an  S-curve  pitching  moment.  The  modified  configuration 
was  assumed  to  have  the  same  lift  and  trim  characteristics  as  the  basic 
cylindrical  afterbody  bomblet,  but  the  subsonic  drag  coefficient  was  in¬ 
creased  by  an  order  of  magnitude. 

Figure  21  compares  the  planar  trajectory  and  range  dispersion  for 
the  basic  and  high-drag  configurations,  based  on  the  representative  subsonic 
dispenser  event  condition.  Bomblet  orientation  for  maximum  dispersion  is 
assumed.  The  trajectory  data  show  that  a  large  amount  of  dispersion  is  lost 
due  to  the  higher  drag,  even  though  the  flight  time  is  increased. 

The  increased  drag  was  effective  in  increasing  the  maximum  flight 
path  angle  at  impact  from  65  to  81  degrees  and  the  minimum  flight  path 
angl  -  at  impact  from  10  degrees  to  72  degrees. 

Dual -Mode  Bomblet  The  dual-mode  S-Curve  bomblet  incorporates 

a  drag  device  which  can  be  deployed  after  a  pre-selected  time*  of  flight. 
Dispersion  is  achieved  during  the  first  (lifting)  flight  phase:  the 

second  (high-drag)  flight  phase  is  utilized  for  increasing  the  flight  path 
angle  at  impact. 

Trajectory  and  dispersion  data  for  the  dual -mode  type  R- Curve 
bomblet  are  compared  with  the  standard  S-Curve  bomblet  in  Figure  22.  Re¬ 
sults  for  the  dual -mode  bomblet  are  shown  for  two  different  values  of  the 
drag  coefficient  for  the  high-drag  phase.  The  smaller  drag  coefficient  is 
representative  of  an  AID  configuration  while  the  larger  drag  coefficient  U 
representative  of  a  small  parachute. 

Deployment  of  a  drag  device  at  about  2.0  seconds  after  event  and  at 
about  1000  (eet  above  ground  level  significantly  increases  the  impact  (light 
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path  angle  but,  at  the  same  time,  cause  a  large  loss  in  range  dispersion, 
particularly  for  the  parachute  device. 

To  achieve  both  large  dispersion  and  large  impact  flight  path  angles 
it  is  evident  that  decelerator  deployment  must  be  delayed.  In  fact,  for  the 
lift-down  trajectories,  it  is  net  even  necessary  for  the  decelerator  to  deploy, 
since  the  impact  flight  path  angles  are  large  in  any  case.  With  the  para¬ 
chute-type  decelerator  deployed  at  4,0  seconds  the  range  dispersion  is  about 
1500  feet  while  the  impact  flight  path  angle  is  of  the  order  of  70  degrees  for 
both  the  lift-up  and  lift-down  cases. 

G.  ANALYSIS  OF  ALTERNATE  S-CURVE  CONFIGURATIONS 

Many  bluff- shaped  bodies  with  low  aspect  ratio  fins  or  subcaliber 
stabilizers  can  be  made  to  have  S-Curvc  moment  characteristics  with  either 
small  or  no  modifications. 

Candidate  S-Curve  configurations  of  both  low-fineness  and  high-fine¬ 
ness  ratio  can  be  selected  with  good  lift-curve  slopes.  However,  the  after¬ 
body  configuration  is  quite  important  and  most  boattail  shapes,  while  having 
good  moment  characteristics,  tend  to  experience  a  reduction  in  lift-curve 
slope. 


A  number  of  recent  wind  tunnel  tests  of  bomb  and  bomblet  models 
provide  data  on  configurations  which  differ  significantly  from  the  two  basic 
S-Curve  bomblets  analyzed  in  the  present  studies.  References  6,  7.  8,  14, 
17,  and  18  present  aerodynamic  data  on  configurations  with  S-Curve  mo¬ 
ment  characteristics,  which  may  be  useful  in  the  development  of  alternate 
S-Curve  configurations.  Normal  force  -slopes  for  a  number  of  candidate 
S-Curve  configurations  are  presented  m  Figure  £3. 

It  is  noteworthy  that  S-Curve  versions  of  the  AJ*AM  and  HLU - 87 
bomblets,  b«th  possess  poor  lift  curve  slopes.'* 

Low  -Fineness-  Ratio  Configurations  Low  fineness  ratio  bomblets 
are  sometimes  of  interest  because  of  warhead  compatibility  and  dispenser¬ 
packaging  considerations.  The  practicability  of  a  low- fineness- ratio 
S-Curve  bomblet  has  been  reviewed  and  some  6 -DOF  motion  simulations 
accomplished.  Two  aspects  of  the  pitching  moment  behavior  of  low* -fineness- 
ratio  shapes  are  detrimental  to  good  Highl  characteristics:  first  the  low 
fineness  ratio  body  has  a  small  moment  curve  slope  at  zero  angle  of  attack 
(or  trim  angles  less  than  about  .1U  degrees;  secondly,  these  configurations 


**  S-Curve  type  dispersion  has  been  attempted  with  both  the  APAM  and 
IU.U-87  i References  6  and  19). 
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tend  to  have  a  rearward  center  of  pressure  shift  at  zero  angle  of  attack  with 
increasing  Mach  number,  which  in  many  instances  maker,  the  configuration 
stable  at  small  angles  of  attack  above  the  critical  Mach  number. 

.  Although  the  dispersion  performance  of  the  low- fineness- ratio 
S-Curve  configuration  at  subcritical  Mach  numbers  appears  promising  under 
ideal  conditions,  the  effects  of  u  lintentional  roll  and  initial  disturbances 
tend  to  greatly  reduce  the  dispersion.  One  of  the  basic  problems  is  that  a 
low-- fineness- ratio  configuration  becomes  gyroscopically  stable  at  small  roll 
rates,  as  a  result  of  both  the  reduction  in  the  polar  to  transverse  moments 
of  inertia  ratio  and  the  reduction  in  the  pitching  moment  derivative  with  de¬ 
creasing  fineness  ratio.  For  a  large  diameter  bomblet,  which  was  investi¬ 
gated  in  detail,  gyroscopic  stability  was  achieved  with  a  roll  rate  as  small 
as  3  cps.  Even  at  roll  rates  less  than  that  required  for  gyroscopic  stabili¬ 
zation,  the  gyroscopic  precession  resulting  from  roll  and  the  out-of-trim 
pitching  moment  causes  spiralling  flight. 

It  is  thus  concluded  that  considerable  effort  would  be  required  to 
develop  a  practicable  low -fineness- ratio  S-Curve  bomblet,  and,  probably, 
severe  restrictions  would  have  to  be  imposed  upon  the  release  environment. 

H.  AIR-GUN  LAUNCHED  MODEL  FREE- FLIGHT  TESTS 

Free- flight  tests  of  the  S-Curve -type  bomblet  were  accomplished  at 
test  area  13-82  at  Eglin  Air  Force  Rase,  Florida  during  October  1972.  A 
svS-ineh-bnre  air  gun  from  the  BARS  test  factlity  was  used  for  launching  the 
models  at  nominal  velocities  of  3*0  to  400  fps  and  at  an  initial  flight  path 
angle  of  30  degrees.  The  first  series  of  tests  was  accomplished  with  3-inch- 
diameter  S-Curve  bomblet  models  of  the  cylindrical  afterbody  configuration 
13 SN SdASl  ^  S2 *  These  models  were  launched  Singly  at  zero  angle  of  attack. 
The  3 -inch -diameter  models  were  divided  into  two  groups;  the  first  group  of 
29  models  was  ballasted  nose  heavy  fur  non-lifting  ballistic  flight,  while  the 
second  group  of  32  models  had  the  center  «»f  gravity  Lr»  calibers  from  the 
nose  (corresponding  to  the  basic  S-Curve  configuration).  Five  models  of  the 
second  group  incorporated  intentional  fin  misalignment.  The  second  series 
of  teat*  utilised  l.i -inch -diameter  models,  which  were  launched  in  clusters 
of  seven  by  the  use  of  a  special  sabot.  The  second  series  of  tests  utilised 
36  cylindrical  afterbody  bomblet  models  and  34  boattail  afterbody  bomblet 
models.  Physical  characteristics  data  for  the  flight-test  models  are  sum* 
mart  zed  in  Table  V. 

Velocity  and  altitude  measurement*  for  most  of  the  tests  were 
obtained  from  high-speed  motion  picture  coverage  of  an  8-foe?  segment  of 
the  India*  trajectory.  For  the  standard  3-meh  models,  the  mean  launch 
velor  tty  was  determined  to  be  3f»9  ips  w;th  a  standard  deviation  of  to  fpv. 

For  the  cluster  rounds,  the  mean  velocity  was  360  fps  with  a  standard  devia¬ 
tion  ef  (ps. 
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TABLE  V.  PHYSICAL  CHARACTERISTICS  DATA  FOR  S-CURVE 
FREE-FLIGHT-TEST  MODELS 


Par  amete  r 

Configurations 

Cylindrical  Afterbody 

i - 

Roattail 

Afterbody 

3  - In  ch 

Diameter 

1 . 5  -  I  n  c  h 

Diameter 

1 . 5 -Inch 

Diameter 

Balli  stic 

Standard 

Standard 

Standard 

Model  length  (Inches) 

1  2.  0 

12.  0 

6.  0 

6.  0 

Weight  (Pounds) 

5.  1  3 

5.  13 

1  .  04 

1  .  04 

CG  (Inches  from  nose) 

3.  0 

4.  48 

2.  3 

2.  79 

CG  (Calibers  from  nose) 

1.  00 

1.  49  3 

1  .  5  3 

1  .  86 

Moments  of  Inertia 

Polar  (Slug-ft^) 

0. 00136 

N/A 

N/A 

Transverse  (Slug -ft1*) 

0.  0105 

0. 00874 

N/A 

N/A 

The  lateral  velocity  perturbations  for  the  standard  3-inch  S-Curve 
models  were  of  the  order  of  9  fps  (standard  deviation).  Wind  velocities  dur¬ 
ing  the  tests  were  generally  less  than  5  mph  and  no  wind  corrections  were 
applied  to  the  test  impact  pattern  data. 

Impact  Patterns  Flight  test  impact  patterns  for  the  3-inch  diam¬ 
eter  models  are  shown  in  Figures  24  and  25.  Figure  24  depicts  the  impact 
pattern  of  the  non-lifting  ballistic  models,  while  Figure  25  illustrates  the 

impact  pattern  of  the  standard  bomblet. 

\ 

The  test  results  arc  compared  with  Monte  Carlo  simulations  using  . 
the  aerodynamic  coefficient  data  of  Table  I  with  appropriate  adjustment  for 
longitudinal  center  of  gravity.  The  Monte  Carlo  simulations  also  incorporate 
the  velocity  perturbations  described  above.  For  the  Monte  t  arlo  simula¬ 
tions  the  configurational  asymmetries  were  assumed  to  be  the  same  as  for 
the  prototype  bomblets  (paragraph  E),  including  a  standard  deviation  value 
for  eg  lateral  offset  of  0.0016  inch. 

The  test  impact  pattern  for  the  ballistic  models  (Figure  24)  shows 
relatively  small  cross- range  dispersion,  as  would  be  expected.  The  mean 
center  of  impact  as  determined  from  the  test  data  is  3110  feet  from  the 
launch  point. 

The  Monte  Carlo  simulation  for  the  ballistic  models  results  in  a 
mean  center  of  impact  at  3310  feet,  or  200  feet  greater  range  than 
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Figure  24,  Flight  Test  ami  Simulated  Impact  Patterns  for 

Atr-Ciun  Launch*-  !  i-inch  Diameter  Ballistic  Models 
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CROSS-  RANGE  (Ft#t) 


RANGE  FROM  GUN  (Feet) 

b)  Flight  Test  Data 

Figure  2r*.  Flight  Teut  and  Simulated  Impact  Patterns  tor 

Air-Gun  Launched  l-inch  Diameter  Standard  S- Curve 
Lombiet  Models  (Cylindrical  Afterbody) 


the  test  data.  The  predicted  cross-range  dispersion  also  exceeds  the  ob¬ 
served  cross-range  dispersion.  Thus,  the  cross-range  velocity  perturba¬ 
tions  used  in  the  simulations  (which  were  derived  from  the  photo  coverage 
of  the  test  launchings)  are  probably  too  large. 

The  test  impact  pattern  for  the  standard  3 -inch  S-Curve  models 
(Figure  25)  confirms  the  effectiveness  of  the  S-Curve  dispersion  concept. 

The  test  impact  pattern  is  approximately  800  feet  in  width  and  1200  feet  in 
length,  with  the  mean  center  of  impact  approximately  2800  feet  from  the 
launch  point. 

The  Monte  Carlo  simulation  for  the  3-inch  S-Curve  models  results 
in  an  impact  pattern  similar  to  the  test  results,  with  the  mean  center  of 
impact  about  3000  feet  from  the  launch  point.  The  standard  deviations  of 
range  and  cross  range,  as  determined  for  both  the  flight-test  and  simulated 
impact  patterns,  are  summarized  in  Table  VI.  Comparison  of  the  test  data 
with  the  Monte  Carlo  simulation  shows  that  the  range  dispersion  is  slightly 
under-predicted  while  the  cross-range  dispersion  is  over-predicted.  The 
initial  cross-range  velocity  perturbations  (which  are  probably  too  large) 
contribute  to  the  over-prediction  of  the  crosi  range. 

The  effect  of  intentional  fin  misalignment  on  the  test  models  is  also 
indicated  in  Figure  25.  Although  the  sample  is  small,  an  improvement  in 
cross -range  dispersion  is  noted  with  intentional  fin  cant  and  roll. 

A  group  of  10  recovered  3-inch-diameter  S-Curve  models  were 
re -flown  to  determine  the  effect  of  nose  scratches  and  other  asymmetries 
resulting  from  impact  damage.  Dispersion  data  for  this  test  group  are  sum¬ 
marized  in  Table  VI.  As  expected,  the  test  data  indicate  a  moderate  de¬ 
crease  in  both  dispersion  and  range. 

Test  results  for  the  1.5 -inch -diameter  cluster- launched  models  are 
shown  in  Figure  26,  An  extremely  large  impact  pattern  was  achieved  for 
the  boattail  bomblet,  while  the  cylindrical  afterbody  bomblet  displayed  an 
elongated  pattern  with  reduced  cross -range  dispersion  compared  with  the 
.1  -inch- diameter  cylindrical -afterbody  bomblet. 

Although  Monte  Carlo  simulations  were  not  attempted  for  the  cluster 
shots  (because  the  initial  conditions  were  poorly  known),  it  was  anticipated 
that  the  dispersion  would  increase  in  proportion  to  the  trimmed  Hft/weigat 
ratio  for  the  respective  model  configurations.  Thus,  the  dispersion  of  the 
1.5-inch -diameter  cylindrical  and  boattail -afterbody  models  was  expected 
to  increase  compared  to  the  3-inch -diameter  cylindrical-afterbody  model 
by  the  lift/weight  ratios  which  were  1,23  and  1.44,  respectively.  The  cross  - 
range  dispersion  of  the  1.5-inch-diameter  boattail  models  is  in  approximate 
agreement  with  the  lift/wcight  ratio,  but  the  cross- range  dispersion  of  the 
1.5- inch -diameter  cylindrical-afterbody  models  does  not  agree  with  the 
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lift/weight  ratio.  No  explanation  is  yet  available  for  the  unusual  perform¬ 
ance  of  the  1.5-inch-diameter  cylindrical-afterbody  models. 


TABLE  VI.  SUMMARY  OF  IMPACT  DISPERSION  DATA  FOR 

AIR -GUN -LAUNCHED  FREE- FLIGHT  MODELS 

^ean 

Range 

(Feet) 

Standard 

Deviation 

Range 

(Feet) 

Standard 

Deviation 

Cross-Range 

(Feet) 

3  -  Inch-  Diameter  Ballistic  Models 

t 

Test  Results 

3110 

271 

27 

Monte  Carlo  Simulation 

3310 

105 

119 

3-Inch- Diameter  Standard 

S-Curve  Models 

Test  Results  a 

2783 

249 

177 

Monte  Carlo  Simulation 

3033 

222 

322 

Test  Results  for  Re- Flown 
Models 

2542 

227 

122 

1,5- Inch -Diameter  S-Curve  Models 

Test  Results- -Cylindrical 
Afterbody 

1981 

333 

108 

Test  Results --Boattail 

Afterbody 

2033 

642 

233 

*  Includes  models  with  intentional  fin  misalignment. 


SECTION  m 


ROLL- THROUGH- ZERO  BOMBLET  INVESTIGATIONS 


A,  CONFIGURATIONS 


During  the  present  contractual  effort  attention  was  focused  on  the 
development  of  a  practical  roll-through-zero  bomblet  configuration.  The 
principal  configurational  requirements  for  this  type  bomblet  are  (1)  a  large 
trimmed  lift  coefficient,  (2)  a  means  for  trimming  the  bomblet  in  pitch,  and 
(3)  a  means  for  reversing  the  initial  direction  of  roll. 


The  roll  reversing  torque  must  be  small,  since  the  deflection  of  the 
bomblet  trajectory,  and  hence  the  dispersion,  has  previously  been  shown 
(Reference  1)  to  be  related  to  the  static  roll  torque  coefficient  by  the  expres 
sion 


AY 


CN 

m' 


(2) 


where  AY  is  the  total  deflection  of  a  roll-through-zero  trajectory  in  rela¬ 
tion  to  a  corresponding  ballistic  trajectory. 


In  addition,  the  bomblet  must  have  a  relatively  large  stability  margin 
at  trim  such  that  the  pitch  natural  frequency  is  above  the  intended  range  of 
roll  rates.  This,  in  turn,  requires  that  the  trim  moment  be  relatively  large. 
At  the  same  time,  the  induced  roll  moments  due  to  the  trim  surfaces  must  be 
small  in  comparison  with  roll  reversing  torque  coefficient,  .  When  all 

of  these  factors  are  considered,  the  design  of  the  roll-through-zero  bomblet 
becomes  quite  complex. 

Investigation  has  shown  that  wings  and  canards,  while  offering  good 
trimmed  lift,  tend  to  have  highly  unfavorable  induced  roll  and  side  moments, 
precluding  their  use.  Likewise,  estimates  have  shown  that  a  d-ag-type  trim 
device  will,  in  general,  provide  insufficient  trim  moment,  and  may  also  lead 
to  sizable  induced  roll  and  side  moments  through  interference  effects. 


The  use  of  intentional  fiq  incidence  for  trim  was*  proposed  as  the  re¬ 
sult  of  preliminary  studies  (Reference  l).  During  the  present  program  the 
aerodynamics,  flight  dynamics,  and  dispersion  characteristics  of  this  type 
configuration  were  extensively  investigated.  To  allow  direct  comparison  of 
the  performance  of  the  roll -through -zero  bomblet  with  the  S-Curve-type 
bomblet,  the  body  and  i>n  geometry  were  made  identical  to  the  S- Curve  con¬ 
figuration  .  f  in  r.i-dence  of  10  degrees  was  then  applied  to 

two  opposing  fins,  and  an  effective  fin  cant  angle  of  0.1  degree  was 


•>0 


introduced.  The  resulting  basic  roll-through-zero  configuration  is  depicted 
in  Figure  27. 

The  roll-through-zero  bomblet  configuration  discussed  in  this  report 
is  scaled  to  a  diameter  of  two  inches;  and  the  values  for  the  weight  and 
inertia  are  assumed  to  be  identical  to  those  for  the  S-Curve-type  bomblet, 
namely, 

W  =  1.51  lbs 

Ix  ^  1.629  x  10‘4  slug-ft2 

I  =  1.06  x  10“ ^  slug-ft2 

In  general,  it  is  assumed  that  the  roll -through- zero  bomblet  will 
have  canted  fins  and  will  be  released  from  a  spinning  dispenser  with  a  known 
roll  rate  of  the  order  of  1  to  2  cps.  However  the  requirement  for  a  specific 
initial  roll  rate  and  roll  direction  is  not  mandatory  since  it  is  possible  for 
the  bomblet  to  provide  its  own  roll  rate  reversal.  Figure  28  shows  sche¬ 
matically  a  roll  reversing  fin  tab,  which  functions  with  decreasing  dynamic 
pressure.  Other  devices,  which  function  on  the  basis  of  a  delay  time,  are 
easily  envisioned. 

B.  AERODYNAMIC  CHARACTERISTICS 

Static  Aerodynamics  The  static  aerodynamic  characteristics  of  the 
basic  roll-through- zero  bomblet  configuration  were  evaluated  through  the 
transonic  range  in  the  4T  wind  tunnel  at  AEDC.  For  each  Mach  number  and 
angle  of  attack  six  component  force  and  moment  data  were  measured  for 
aerodynamic  roll  angles  from  -165  to  180  degrees  at  intervals  of  15  de¬ 
grees.  The  trim  roll  orientation  for  positive  lift  corresponds  to  an  aero¬ 
dynamic  roll  angle  of  180  degrees.  Aerodynamic  coefficient  plots  versus 
aerodynamic  roll  angle  for  angles  of  attack  of  10,  15,  and  20  degrees  and 
Mach  number  0.5  are  shown  in  Figure  29. 

The  large  effect  of  the  aerodynamic  roil  angle  is  clearly  evident 
when  the  coefficients  are  plotted  as  a  function  of  $.  The  first  harmonic  in 
roll  with  period  of  2  *  clearly  predominates  and  is  due  to  the  fin  incidence. 
Surprisingly,  the  second  and  third  harmonic  arc  also  evident,  with  periods 
of  if  and  2 it/ 3,  respectively,  while  the  usual  sixth  harmonic  (due  to  the 
presence  of  six  fins)  is  not  apparent. 

The  data  of  Figure  29  also  show  the  loss  in  normal  force  due  to  the 
trim  incidence,  and  the  unstable  roil  moment  variation,  d  Cj/d4»  >  0,  at 
the  trim  orientation  is  noteworthy. 

For  bomblet  design  purposes,  it  is  important  to  have  a  quantitative 
measure  of  the  fin  normal  force  effectiveness  due  to  incidence.  The  ratio 
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Figure  11.  Basie  RoU«  Through -/.tiro  Bomblet 
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Figure  28.  Roll  Reversing  Tab 


/Cjvj  was  determined  both  from  the  experimental  data  and  estimated 

vising  the  method  of  Reference  20  in  conjunction  with  fin  aerodynamic  force 
data  from  Reference  21.  The  theory  was  found  to  underpredict  the 

/CN  ratio  by  a  factor  of  about  two  at  zero  angle  of  attack,  but  at  a 

10-degree  angle  of  attack  the  experimental  and  theoretical  values  for 

were  0.58  and  0.53,  respectively.  The  highly  nonlinear  varia¬ 
tion  of  fin  effectiveness  with  angle  of  attack  complicates  the  theoretical 
aerodynamic  design,  and  wind  tunnel  evaluation  of  the  bomblet  trim  charac¬ 
teristics  are  nearly  essential. 

To  obtain  a  nominal  trim  angle  of  attack  of  10  degrees  at  low  sub¬ 
sonic  velocities,  the  bomblet  longitudinal  center -of- gravity  was  positioned 
1.185  calibers  from  the  nose.  A  more  rearward  eg  would  be  possible  with 
increased  span  fins,  and  the  aerodynamic  characteristics  of  such  a  configu¬ 
ration  can  be  determined  from  the  test  data  of  Reference  13. 

Dynamic  Stability  Parameters  Pitch  damping,  magnus,  fin  roll 
effectiveness,  and  roll  damping  data  for  the  roll-through-zero  bomblet  are 
based  on  test  results  for  the  S-Curve-  bomblet  configuration  BgNg2AgjFg2» 

Aerodynamic  Coefficients  for  6 -DOF  Trajectory  Program  Adapta¬ 
tion  of  the  wind  tunnel  aerodynamic  coefficients  to  the  6 -DOF  trajectory 
program  (Reference  16)  was  accomplished  using  first  and  second  harmonics 
for  Cjyj,  Cj^,  and  Cgp  and  first,  second  and  third  harmonics  for  Cgj^  and 
Cj^  .  The  first  harmonics  were  introduced  through  the  coefficients  Cjsj  ^ 

and  which  are  related  to  Cjsj,  Cgp,  C^,  and  Cgj^  through  the  aero¬ 

dynamic  roll  angle,  while  higher  harmonics  in  roll  were  incorporated  by  the 
coefficients  C^jj,  CMj,  Cgpj,  CgM2,  Cgj^,  C C  as  de¬ 

fined  in  Reference  16.  A  typical  fit  to  the  experimental  coefficients  using 
th  ree  harmonics  is  illustrated  in  Figure  30. 

Table  VII  summarizes  the  aerodynamic  coefficients  used  in  the 
6 -DOF  motion  studies  and  Monte  Carlo  impact  pattern  simulations. 


C.  PREDICTED  FLIGHT  DYNAMICS 

Figure  31  illustrates  a  typical  set  of  motion  histories  for  the  basic 
roll -through* zero  bomblet  configuration,  based  on  the  standard  subsonic 
event  conditions  where  HOB  =  2000  feet,  -  900  fps,  and  Y£  -  45  degrees 
The  data  are  selected  from  one  of  a  series  of  trajectories  comprising  a 
Monte  Carlo  dispersion  study.  This  particular  bomblet  has  an  initial  nega¬ 
tive  roll  rate  of  one  rps,  a  mean  effectiveness  fin  cant  of  0,115  degree,  and 
a  eg  lateral  offset  of  0.001  inch,  it  is  to  be  noted  that  for  this  particular 
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TABLE  VII.  AERODYNAMIC  COEFFICIENT  SUMMARY  FOR  BASIC 
ROLL- THROUGH-  ZERO  BOMB  LET  CONFIGURATION 


bomblet  the  maximum  induced  roll  moment  due  fo  eg  lateral  offset  is  about 
60  percent  of  that  due  to  the  intentional  fin  cant. 

The  variables  used  to  describe  the  motion  are  the  angle  of  attack 
plane  orientation,  $  ,  the  aerodynamic  roll  angle,  <*>  ,  the  roll  rate,  p  ,  and 
the  total  angle  of  attack,  a  . 

Examination  of  the  motion  histories  shows  that  about  0.5  second  are 
required  for  the  bomblet  to  attain  zero  roll  rate  and  for  the  roll  angle  to 
stabilize  at  the  trim  orientation.  The  initial  motion  is  well- behaved,  with 
the  angle  of  attack  plane  remaining  in  a  single  quadrant  for  about  one  second 
of  flight  such  that  large  dispersion  can  occur. 

During  the  initial  phase  of  flight  and  until  about  three  seconds  after 
release  the  aerodynamic  roll  angle  remains  stable  at  =  it  ,  even  though 
the  total  aerodynamic  induced  roll  torque  has  an  unstable,  or  positive, 
gradient  with  respect  to  the  aerodynamic  roll  angle.  The  stable  motion  be¬ 
havior  can  be  explained  by  the  side  moment,  which  at  the  trim  angle  of 
attack  has  a  maximum  value  an  order  of  magnitude  larger  than  the  maximum 
induced  roll  moment.  At  trim  the  gradient  of  the  side  moment  Is  such  that 
the  coning  motion  tends  to  stabilize  the  trim  orientation.  ^ 

After  about  three  seconds  of  flight  the  nature  of  the  motion  changes 
due  to  the  increasing  roll  rate,  which  is  approaching  the  bomblet  pitch-roll 
resonance  frequency.  This  results  in  the  familiar  trim  phase  shift  and  trim 
amplification,  both  of  which  effects  are  easily  seen  in  the  motion  data. 

The  necessity  for  suitably  restricting  the  eg  lateral  offset  of  the 
roll -through- zero  bomblet  must  be  emphasized.  Because  the  aerodynamic 
roll  angle  undei  trim  conditions  is  invariant  with  respect  to  the  incidence 
plane,  random  orientation  of  the  eg  lateral  offset  will  produce  roll  torques 
which  are  both  random  in  magnitude  and  direction,  with  attendant  large  var¬ 
iation,'.  in  roll  torque  between  bomblet s.  These  roll  torque  variations  could, 
for  even  a  moderately  large  eg  offset,  exceed  the  intentional  fin  cant  roll 
torque,  and  result  in  rapid  roll  rate  butld-up  to  values  exceeding  the  pitch 
resonance  frequency. 

D.  Dl  S  PH  H  SI  ON  P  HE  DICTION 

As  with  the  S-Curve  bomblet,  the  intent  of  the  present  effort  has  been 
to  provide  realistic  impact  dispersion  simulations,  using  the  Monte  Carlo 
method.  To  provide  for  comparison  o*  the  S-Curve  and  roll -through- zero 


4 

This  fact  was  proven  by  reversal  of  the  roll  moment  and  side  moment 
gradients.  With  a  stable  roll  moment  and  unstable  side  moment,  the 
motion  became  violently  unstable  with  no  tendency  for  trim. 


bomblet  performance,  the  standard  subsonic  event  condition  was  utilized  as 
the  basis  for  the  dispersion  simulations. 

The  physical  characteristics  of  the  2-inch-diameter  roll-through- 
zero  bomblet  have  been  previously  described.  For  the  dispersion  simula¬ 
tions,  the  mean  fin  cant  angle  was  assumed  to  be  0.100  degree  with  a 
standard  deviation  m  error  of  0.0816  degree.  The  standard  deviation  of  eg 
offset  was  taken  to  be  0.0016  inch.  Because  of  the  large  fin  incidence  em¬ 
ployed  for  trim,  no  other  lateral  misalignments  were  considered. 

The  computed  impact  pattern  for  the  roll-through- zero  bomblet  is 
shown  in  Figure  32.  The  pattern  is  of  approximately  circular  shape  with  a 
diameter  of  about  hhO  feet.  The  standard  deviations  of  range  and  cross- 
range:  dispersion  are  summarized  in  Table  VIII,  along  with  the  comparative 

data  for  the  S- Curve  bomblet. 


TABLE  VIII.  COMPARISON  OF  IMPACT  DISPERSION  DATA  FOR 
ROLL -THROUGH-ZERO  AND  S-CURVE  BOMBLETS  - 
SUBSONIC  EVENT  CONDITIONS 


Homblot  f yp‘‘ 

Mean 

Range 

X 

( Feet) 

Mean 

Cross-Range 

Y 

(Feet) 

Standard 

Deviation 

Range 

°X 

( Feet) 

Standard 

Deviation 
Cross-  Range 

0  Y 

( Feet) 

l<>. l!  -  Through 

t  805 

-28 

1  33 

1  59 

S  Curve 

1899 

-26 

246 

185 

Cornua  ns  cm  with  S-Curve  Bomblet  The  roll-through-zero  bomblet 
configuration  considered  in  this  investigation  has  effective  dispersion,  but 
the  area  coverage  is  not  quite  as  large  as  that  achieved  by  the  5- Curve  type 
in, mulct.  The  primary  reason  for  the  smaller  area  covera-  ■  ’  the  roll  - 

through-  zero  bomblet  is  the  lower  value  of  trimmed  lift,  ALL-  tgh  both 
types  uf  buinbiet*  have  been  designed  for  a  nominal  trim  angle  of  attack  of 
at, out  in  tu-groes.  the  lift  curve  slope  of  the  present  roll- through -zero 
bomblet  is  appreciably  reduced  by  the  fin  negative  lift  force  required  for 
mm.  At  Mach  number  0.5  and  at  10  degrees  angle  ot  attack,  the  roll- 
through  -  r.f-rc;  bomblet  has  a  lift  coefficient  of  0.45,  compared  to  a  lift  co- 
t-fficirn!  C.5S  for  the  S-Curve  bomblet.  With  allowance  for  the  lift 
coefficient,  the  dispersion  of  the  roll  -  through  -  zero  and  £  --Curve  bumble!* 

,«  eomoarable.  It  *  significant  that  the  trimmed  lift  coefficient  for  the 
ro?5- through- /.rr<»  *-•  mblet  can  be  increased  to  a  value  of  0.#^  at  10  degrees 
4 ui  attack  by  r  -  Teasing  the  fin  exposed  semi -span  from  O.CVd  to  0. 140 
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Figure  32.  Monte  Carlo  impact  Pattern  Prediction  for  Basic 

Roll  -  Through- /.ero  Bomblet  -  High-Subsonic  Event 
Condition 


Comparison  with  Theory  A  simplified  theoretical  prediction  of  the 
roll -through- *ero  bomblet  dispersion  can  be  obtained  using  equation  (2)  in 
conjunction  with  the  dispersion  approximations 


r 


h 

sin  > 


AY 


(3) 


(4) 


where  r  and  R  represent  thr  >  ro.s?.. range  dispersion  and  range  dispersion, 
respectively,  at  ground  impact,  in  the  theoretical  calculations,  it  is  as¬ 
sumed  that  the  roll  torque  coefficient.  C£() ,  is  due  only  to  the  mean  fin  cant 
angle.  The  induced  roll  torques  are  not  included  because  they  are  random. 
Table  IX  compares  the  standard  deviation  values  of  range  and  cross- range* 
dispersion,  as  obtained  from  the  Monte  Carlo  simulations,  with  the  theoret¬ 
ical  predictions  for  r  and  R. 
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TABLE  IX.  COMPARISON  OF  THEORETICAL  AND  MONTE  CARLO 
DISPERSION  PREDICTIONS  FOR  ROLL- THROUGH- ZERO 

BOMBLET 

_ 

Cross-Range 

Di  spersion 
(Feet) 

Range 

Dispersion 

(Feet) 

Monte  Carlo;  o„  and  o 

X  y 

1  59 

133 

Theoretical;  r  and  R 

540 

764 

The  reason  for  the  poor  agreement  between  the  theoretical  and 
Monte  Carlo  predictions  is  believed  to  be  due  primarily  to  the  fact  that  the 
induced  roll  torques  are  not  included  in  the  theoretical  prediction. 
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SECTION  IV 


DEVELOPMENT  OF  MONTE  CARLO  METHOD 
FOR  DISPERSION  PREDICTION 

A.  GENERAL 

The  Monte  Carlo  technique  is  a  procedure  by  which  one  can  obtain 
approximate  evaluations  of  mathematical  relationships  which  involve  one  or 
more  probability  distribution  functions. 

The  Monte  Carlo  technique  consists  of  simulating  an  experiment  to 
determine  some  probabilistic  property  of  a  population  of  events  by  the  use 
of  random  sampling  applied  to  the  components  of  the  events.  In  the  present 
application  of  this  mathematical  definition,  the  experiment  is  the  deploy¬ 
ment  of  a  bomblet  cluster,  the  events  are  the  bomblet  trajectories,  a  prob¬ 
abilistic  property  of  the  events  is  the  impact  pattern,  and  the  random  sam¬ 
pling  is  applied  to  those  variables  which  affect  the  bomblet  trajectories, 
namely,  the  initial  motion  parameters  and  the  bomblet  configurational 
asymmetries.  We  could  also  include  atmospheric  or  flow  disturbances 
along  the  trajectories. 

The  application  of  the  Monte  Carlo  technique  to  the  flight  period 
immediately  after  dispenser  opening  is  extremely  difficult  because  during 
this  period  multi -component  flow  conditions  and  very  large  bomblet- 
bomblet  and  bomblet-di spenser  interference  effects  exist,  and  these  have 
to  be  describable  in  probabilistic  terms  for  the  trajectories  to  be  simulated. 
However,  after  a  brief  interval  time,  break-up  of  the  bomblet  cluster 
usually  occurs,  such  that  some  of  the  bomblets  are  free  of  interference 
effects.  If  the  probable  flight  conditions  at  cluster  break-out  are  estab¬ 
lished,  then  the  bomblet  trajectories  can  be  initiated  at  the  time  of  break¬ 
out,  and  the  trajectory  computations  are  greatly  simplified. 

B.  MONTE  CARLO  TRAJECTORY  PROGRAM 

Using  the  cluster  break-out  concept,  a  Monte  Carlo  trajectory  pro¬ 
gram  was  developed.  This  program  provides  for  the  sequential  calculation 
of  a  pre-set  number  of  6-DOF  trajectories  and  corresponding  impact  points, 
using  statistical  input  for  the  cluster  break-out  conditions  and  configura¬ 
tional  asymmetries.  To  provide  the  required  generality  in  the  equations  of 
motion,  a  number  of  modifications  and  additions  were  made  to  the  Extended 
Capability  Magnus  Rotor  and  Ballistic  Body  6-DOF  Trajectory  Program 
(Reference  31,  and  a  new  subroutine  MONCAR  was  incorporated  to  provide 
statistical  input  data  to  the  6-DOF  trajectory  program.  A  complete  descrip¬ 
tion  of  the  new  Monte  Carlo  trajectory  program  is  provided  in  Reference  16. 
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initial  Motion  Parameters  The  Monte  Carlo  trajectory  program  is 
structured  such  that  the  bomblet  £-DOF  motion  prediction  commences  at 
clu.-ter  break-out.  The  general  assumptions  which  are  used  to  establish  the 
conditions  at  cluster  break-out  are  as  follows: 


bomblet  has  a  time  of  cluster  break-out,  t^  ,  defined  by  an 
unspecitied  probability  distribution. 

2.  Prior  to  cluster  break-out,  the  leading  edge  of  the  cluster 
one  ole  rates  uniformly  and  follows  a  straight  flight  path,  coin¬ 
cident  with  the  dispenser  flight  path  at  event. 

3.  The  cluster  is  radially  symmetric  with  respect  to  bomblet  posi¬ 
tion  and  velocity. 


4.  Cluster  break-out  begins  from  the  cluster  leading  edge  and  on 
the  'duster  axis  of  rotational  symmetry. 


Figure  33  shows,  schematically,  the  variables  which  control  the 
bomblet  initial  position,  velocity,  and  pitch  attitude  at  cluster  break-out. 
The  probabilistic  initial  motion  parameters  and  their  probability  distribu¬ 
tions  lo -  ..pr rifled  event  conditions,  are  given  below. 


Variable 

tu 

v-< 

«  R 


V 


C 


? 


Description 

Cluster  break-out  time 

■tadial  velocity  at  break-out 

Velocity  orientation  at  break-out 

Angle  *»f  attack  at  cluster  break-out 

Angle  of  attack  plane  orientation  at 
cluster  break-out 

Pitch  rate  in  a  plane  at  cluster 
break-out 

Yaw  rate  at  cluster  break-out 

*»ol!  rate  at  cluster  break-out 

Hull  orientation  of  body  axes  at 
cluster  break-out. 


Pro  bability 
Distribution 

Arbitrary 

Arbitrary 

Uniformly  random 

Arbitrary 

Uniformly  random 

Arhlt  rary 

Arbitrary 
Arbitrary 
Uniformly  random 


•  rum  these  quantities,  the  event  conditions  and  the  specified  cluster 
ration.  I1'-  initial  conditions  required  for  the  basic  6-DOF  trajectory 
?,-  ■  df'lcrmtued. 
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VERTICAL 


Figure  Schematic  of  l\omblt*t  Cluster  Hr«?ak-Up  And  Definition 
of  Initial  Motion  Perturbation » 


SrfegiFl'i^wajaMii' 


Configurational  Asymmetries  The  Monte  Carlo  subroutines  and 
6 -DOF  equations  of  motion  have  been  formulated  in  such  a  manner  that  the 
following  geometric  and  inertial  asymmetries  can  be  input  in  probabilistic 


fo  rm, 

Probability 

Symbol 

Description 

Distribution 

A 

Fin  cant 

Normal 

CANT 

X  I, AT 

Lateral  misalignment 

Normal 

Ay 

Lateral  eg  offset 

Normal 

Ax 

Axial  eg  offset 

Normal 

Am 

Mass 

Normal 

tlx 

Axial  moment  of  inertia 

Normal 

A  I 

Transverse  moment  of  inertia 

Normal 

Alxy 

Product  of  inertia 

Normal 

Complete  generality  in  angular  orientation  of  the  asymmetries  is 
assumed  by  introducing  the  following  probabilistic  angles.  ^ 


Symbol 

i 


Probability 

Description  Distribution 

Roll  orientation  of  aerodynamic  surfaces  Uniform 

Roll  orientation  of  aerodynamic  misalignment  Uniform 


Aerodynamic  Asymmetries  Aerodynamic  asymmetry  coefficients 
Cyf).  Cz<>,  Cn,u.  and  Cr,  can  also  be  introduced  directly  into  the  Monte 
Carlo  simulation  but.  in  most  instances,  fin  cant  asymmetry  and  lateral 
misalignment  adequately  describe  both  the  geometric  and  aerodynamic 
asymmetries.  However,  for  some  configurations  it  is  also  necessary  to 
account  for  flow  asymmetries  which  cannot  be  related  to  fin  cant  or  lateral 
misalignment.  Such  flow  asymmetries  can  be  triggered,  for  m.  tance,  by 
almost  undetectable  nose  roughness,  and  can  result  in  so-called  psuedo 
mag  mis  t*»rces  and  moments  as  well  as  induced  roll  moments-  Provision 
ha*  h***.»-  made  in  the  Monte  Carlo  trajectory  program  for  such  phenomena. 


Probability  Distribution  Functions  The  Monte  Carlo  compute1' 
program  provides  for  three  types  of  probability  distribution  functions. 


‘1  be  unbalance,  \  ! 


b'  .»*  •? 


baJ  .»*  <  e 


xy  * 


restricted  to  the  same  plane  as  the 


1.  Uniformly  random 

2.  Normal  distribution 

3.  Arbitrary  distribution  function 

All  of  the  probabilistic  variables  are  determined  from  a  random  number 
subroutine  RANF. 


The  uniformly  random  variables  are  determined  directly  by  multipli¬ 
cation  of  RANF  and  an  appropriate  constant. 


The  normal  distribution  is  obtained  from  the  specified  mean  value  of 
the  variable  and  normal  deviates.  The  normal  deviates  are  obtained  using 
the  Gaussian  approximation 


GAUSS 


-6  i 


12 

RANF  (Xi ) 

1 


(5) 


The  arbitrary  distribution  functions  are  provided  by  tables  where  the 
arguments  are  the  cumulative  frequency  normalized  0  »  1.0  and  the  ordi¬ 
nates  are  the  probable  values  of  the  variable. 


C.  STATISTICAL  ANALYSIS  OF  CLUSTER  BREAK-UP 

To  obtain  the  probability  distribution  functions  for  the  motion  param¬ 
eters  at  cluster  break-out,  it  was  necessary  to  examine  experimental  cluster 
break-up  data.  Photographic  data  from  sled  tests  of  a  Lance  missile  bomb- 
let  warhead  with  approximately  1600  XM-41  bomblets  were  analyzed  to 
determine  cluster  leading  edge  velocity  and  lateral  and  longitudinal  dispersal 
as  a  function  of  time. 

Cluster  Break-Out  Time  Figure  34  illustrates  the  cumulative  fre¬ 
quency  distribution  as  a  function  of  time  of  those  bomblets  which  have  at¬ 
tained  separation  from  the  dense  central  duster  (separation  is  defined  as  a 
density  of  bomblets  of  35  units  per  square  yard  or  less  in  a  sectional  view 
of  the  cluster,  or  about  10  diameters  bomblet- bomblet  separation),  for  a 
warhead  event  velocity  of  19?5  fps.  A  similar  statistical  variation,  adjusted 
to  the  available  data  for  a  subsonic  APAM  bomblet  cluster,  is  also  pre¬ 
sented  in  Figure  34. 

Cluster  Deceleration  Figure  35  illustrates  the  longitudinal  decel¬ 
eration  of  the  duster  leading  edge  as  a  function  of  event  velocity,  based  on 
two  separate  sled  tests  of  the  Lance  warhead.  For  each  of  the  tests,  the 
deceleration  was  found  to  be  nearly  constant  during  the  cluster  break-up 
period.  When  plotted  as  a  function  of  velocity  the  results  show  that  the 


cluster  deceleration  is  very  nearly  proportional  to  the  square  of  the  event 
velocity. 

Radial  Velocity  The  lateral  expansion  of  bomblet  clusters  has  been 
investigated  both  from  experimental  data  and  theoretical  considerations  as  a 
means  of  determining  the  bomblet  radial  velocity  statistics. 

The  Lance  missile  sled  test  data  indicate  maximum  aerodynamic- 
induced  lateral  cluster- spread  velocities  of  the  order  of  130  fps  for  an  event 
velocity  of  2175  fps.  The  frequency  distribution  is  approximately  parabolic 
with  respect  to  a  zero  mean. 

Ground  impact  pattern  data  for  a  non- selfdisper sing  cluster  muni¬ 
tion  (Reference  22)  were  also  analyzed  as  a  means  of  obtaining  the  initial 
radial  velocity  perturbations.  Theoretical  dispersion  patterns  were  com  - 
puted  as  a  function  of  the  initial  radial  velocity  and  matched  to  the  observed 
cross-range  impact  dispersion  to  obtain  an  estimate  of  the  actual  radial 
velocity.  ^  In  a  similar  manner,  the  difference  between  the  predicted  and 
actual  along- range  dispersion  was  used  to  obtain  an  estimate  of  the  cluster 
break-up  time. 

The  above  determinations  of  the  induced  lateral  velocity  are  com¬ 
pared  in  Figure  36  with  measurements  for  APAM  bomblet  clusters  (Refer¬ 
ence  23).  There  is  seen  to  be  very  good  agreement  in  the  two  sources  of 
subsonic  data. 

A  theoretical  prediction  has  been  made  of  the  lateral  velocity  im¬ 
parted  to  a  bomblet  initially  at  the  leading  edge  of  the  cluster.  The  detailed 
development  of  the  theory  is  shown  in  Appendix  1.  The  theory  shows  that  the 
radial  velocity.  Vp  ,  is  a  linear  function  of  the  event  velocity  in  accordance 
with  the  following  formula: 

[R  *1  ^ 

ofrep;  j  <6> 

where  Vg  =  event  velocity 

g  e  acceleration  due  to  gravity 
«-■  s  air  density 


f  It  s.v  noteworthy  that  the  standard  deviation  of  radial  velocity  does  not  cor- 
respond  to  the  standard  deviation  of  cross  range  because  the  mitial  cross¬ 
range  velocity  is  proportional  to  the  sine  of  $■  g,  Assuming  that  is 

distributed,  the  standard  deviation  of  radial  velocity  corresponds 
nr  cumulative  frequency  for  the  cross  range. 


:"-n 
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R  =  cluster  core  rad-. 

(W/CqS)b  =  bomblet  ballistic  ■'efficient. 

If  the  cluster  core  diameter  is  taken  to  be  about  four  times  the  dis¬ 
penser  diameter,  the  theoretical  Vr  is  in  nearly  exact  agreement  with  the 
experimental  data  of  Figure  36.  The  importance  of  the  above  result  is  that 
it  shows  that  the  spread  velocity  is  dependent  upon  both  the  size  of  the 
cluster  and  the  bomblet  ballistic  coefficient. 

Effect  of  Spinning  Dispenser  All  of  the  previous  data  and  analyses 
of  lateral  spreading  are  based  on  small  or  negligible  dispenser  spin  rates. 
For  spinning  dispensers,  an  additional  lateral  velocity  component  will  be 
generated,  which  will  combine  with  the  aerodynamic  induced  velocity,  and 
the  maximum  lateral  velocity  will  be  given  by 

T  1 1/2 

(Vr)totAL  ■  |<VR>AERO  *  <r“>  ]  <7> 

where 

u)  r.  dispenser  angular  velocity 

r  ^  maximum  radial  displacement  of  bomblets 
from  dispenser  axis  of  rotation 

Angle  of  Attack  and  Angular  Velocity  at  Cluster  Break-Out  The 
initial  angle  of  attack  and  angular  velocity  at  cluster  break-up  cannot  be 
generalized,  and  depend  largely  upon  the  bomblet  configuration.  For  small, 
low- fineness- ratio  bodies  without  stabilizing  fins  it  has  been  shown  that  the 
initial  attitude  will  tend  to  be  uniformly  random.  Qualitative  inspection  of 
Roekeye  and  APAM  bomblet  motions  during  cluster  break-up  reveals  that 
the  maximum  angle  of  attack  seldom  exceeds  about  45  degrees.  If  the 
pitching  moment  coefficient  is  known  as  a  function  of  angle  of  attack,  it  is 
possible  to  relate  the  maximum  angle  of  attack  and  the  pitch  rate  at  zero 
angle  of  attack.  It  is  then  reasonable  to  assume  that  the  pitch  rate  at  zero 
angle  of  attack  follows  a  normal  distribution. 

Using  this  procedure  and  assuming  that  the  maximum  angle  of 
attack  (two  sigma)  is  45  degrees,  the  corresponding  pitch  rate  for  the 
S-Curve  and  roll -through -zero  bomblets  is  about  70  radians/second  for 
an  event  velocity  of  000  fps. 


D.  STATISTICAL  ANALYSIS  OF  CONFIGURATIONAL  ASYMMETRIES 

Fin  Cant  Error  The  effective  fin  cant  asymmetry  for  a  multi* 
finned  configuration  can  be  computed  from  alignment  and  accuracy  measure¬ 
ments  on  individual  fins,  or  extracted  from  observations  of  the  spin  behavior 
of  a  representative  quantity  of  test  configurations.  The  resulting  effective 
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fin  cant  statistics  must  be  in  a  form  compatible  with  the  Monte  Carlo  tra¬ 
jectory  program,  which  means  that  the  statistical  variation  of  the  total  fin 
roll  torque  is  required.  The  total  torque,  in  turn,  is  assumed  to  be  pro¬ 
portional  to  an  equal  deflection  of  all  the  fins,  and  this  deflection  is  defined 
as  the  effective  fin  cant. 

By  the  method  of  super  position,  the  total  rolling  moment  due  to 
separate  deflection  of  n  fins  can  be  expressed  as 


Co  =  Co  =  - 

*T  *  TOTAL  n 


(8) 


where  Q,^  is  the  effective  roll  torque  coefficient  due  to  fin  cant,  correspond¬ 
ing  to  simultaneous  deflection  of  all  n  fins  and  6  is  the  effective  fin  cant 
deflection.  Letting  the  variance  of  angular  deflection  of  each  fin  equal  0  p4', 
then  the  standard  deviation  of  C  o_  is 


std  dev 


(9) 


If  the  variance  of  all  fins  is  identical. 


std  dev  C 


£T 


'*6 


_  1/2 


n  o. 


Co  a  v 


*  C*6  (  Vo 


or  (  f  ) 
e  0 


v/V^ 


(10) 


Thu?,  the  effective  fin  cant  angular  error  is  neither  the  standard  deviation 
of  error  for  a  single  fin  nor  an  average  of  the  fin  cant  deviations  of  several 
fins,  but  rather  the  standard  deviation  of  one  fin  divided  by  the  square  root 
of  the  total  number  of  fins. 


To  obtain  representative  values  for  Op,  fin  cant  misalignment  data 
for  81mm  mortor  shells,  2.75-Inch  FFAR  rockets,  Rockeye  bomblets,  and 
Australian  S-Curve  test  bomblets  were  obtained  from  References  24®,  25, 
and  26.  Single  fin  misalignments  for  these  configurations  have  been 


®  The  data  of  Reference  24 
normal  distribution. 


show  that  fin  misalignments  closely  follow  a 
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determined  to  be,  respectively,  0.22  degree,  0.24  to  0,29  degree,  0.2 
degree,  and  0.08  to  0.2  degree.  When  a  range  of  values  is  indicated,  more 
than  one  type  of  configuration  was  inspected.  On  the  basis  of  these  results, 
a  single  fin  alignment  accuracy  of  0.2  degree  appears  to  be  very  reasonable 
for  Monte  Carlo  simulations. 

roi  a  six  fin  S-Curve  or  roil-through-zero  configuration  the  effective 
fin  cant  misalignment  is  0.0816  degree. 

Lateral  Misalignment  Lateral  misalignment  results  primarily 
from  two  sources,  (1)  the  angular  accuracy  of  the  fins  with  respect  to  the 
body  at  their  point  of  attachment,  and  (2)  the  alignment  of  the  fin  assembly 
with  the  remainder  of  the  body.  These  two  types  of  misalignment  tend  to  be 
statistically  independent.  Lateral  misalignment  results  in  transverse  body- 
fixed  forces  and  moments,  in  contrast  to  fin  can.,  which  results  in  pure  roll 
torque.  If  we  define  the  angle  6  as  the  average  misalignment  of  a  planar 
pair  of  fins,  it  is  easily  shown  that  the  standard  deviation  of  effective  6  for 
the  pair  of  fins  is^ 


std  dev  (  6)eff 


°F 

yr  • 


planar  pair 


For  multiple  fins  the  forces  generated  by  (  6  )g£f  are  not  colinear,  and  thus 

the  effective  5  must  be  based  on  the  random  deflection  of  several  pairs  of 
fins.  10  This  :  s  similar  to  the  random  walk  problem  and  leads  to 


std  dev  [  6  ] 

effJm  pairs 


m  pairs 


The  total  lateral  misalignment  is  the  RSS  of  (6)  and  the  fin 

Cl  I 

assembi y- body  misalignment,  c  , 

r  2  2 1 m 

std  dev  (  6  TOTAL>  =  M  std  dev  6  gff)  +  (std  dev  t;  )  j  (13) 

The  plane  of  the  lateral  misalignment  can  obviously  have  any  orienta¬ 
tion  and  is  thus  uniformly  random. 


0  The  definition  of  6  for  lateral  misalignment  is  consistent  with  the  evalua¬ 
tion  of  Cjs] for  a  planar  pair  of  fins. 

J  Fin-fin  interference  effects  are  neglected  here. 
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To  obtain  numerical  values  for  the  effective  lateral  misalignment 
only  the  misalignment  e  need  be  considered  since  it  has  previously  been 
shown  that  a  good  estimate  for  o  jp  is  0.2  degree.  On  the  basis  of  data 
from  References  25  and  26  a  reasonable  value  for  e  is  also  0.2  degree. 
Thus  the  standard  deviation  of  total  lateral  misalignment  is  estimated  to  be 
0.32  degree  for  S-Curve  and  roll -through- zero  bomblets  with  six  fins. 

j 

Static  and  Dynamic  Unbalance  The  static  unbalance  of  bombs, 
rockets,  and  projectiles  depends  greatly  upon  their  design,  manufacturing 
and  loading.  Static  and  dynamic  balance  data  for  ordnance  are  extremely 
limited  and  available  for  only  a  few  configurations.  A  very  comprehensive 
statistical  analysis  of  the  M437,  175mm'projectile  is  presented  in  Refer¬ 
ence  27.  This,  together  with  some  unpublished  data  on  mortar  shells,  was 
the  extent  of  the  available  data.  For  artillery  projectiles  the  standard 
deviation  of  static  unbalance  is  of  the  order  Ay/d  =  0.0004,  and  the  dy¬ 
namic  unbalance  (inclination  of  principal  axis)  is  of  the  order  of  0.001 
degree.  In  contrast,  the  standard  deviation  of  static  unbalance  of  mortar 
shells  is  reported  to  be  as  large  as  Ay/d  =  0.0050. 

Because  of  apparent  large  variations  in  static  unbalance  between 
various  types  of  ordnance,  a  single  value  for  the  standard  deviation  of  Ay, 
applicable  to  all  S-Curve  and  roll -through- zero  bomblets,  was  not  consid¬ 
ered  appropriate.  Consequently,  in  the  Monte  Carlo  simulations  the  static 
unbalance  parameter  is  varied  parametrically  from  zero  to  as  much  as 
0.04  inches. 
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SECTION  V 


DEVELOPMENT  OF  DYNAMIC  WIND  TUNNEL 
MODEL  SUPPORT  SYSTEM 

A*  u  x  ^  l iwivi  Di^oCRIPTION 

The  dynamic  wind  tunnel  model  system  provides  for  testing  of  quasi  - 
axisymmetric  model  configurations  under  simulated  coning  motion  conditions 
at  constant  angle  of  attack,  either  in  steady  or  variable  coning  rate  modes. 
Simultaneously,  the  model  can  be  free  to  roll  about  its  longitudinal  axis  or 
locked-in  at  a  specific  roll  orientation.  The  system  does  not  employ  a  force 
balance;  instead,  the  aerodynamic  moments  are  derived  from  analysis  of 
the  coning  and  rolling  angular  rate  data. 

Design  details  of  the  dynamic  model  support  system  are  presented  in 
Figures  37  and  38.  The  principal  structural  member  is  the  rotating  sting, 
which  is  supported  by  low  friction  bearings.  The  offset  portion  of  the  sting 
can  be  laterally  adjusted  by  means  of  a  boring  head  adapter.  Provision  is 
also  made  for  attachment  of  static  and  dynamic  balance  weights  to  the  shaft 
at  two  stations.  The  sting  is  connected  to  a  Gast  type-4  AM  air  motor  with 
an  output  of  0.8  HP  at  1000  rpm  with  100  psi  supply  pressure.  The  sting 
and  air  motor  are  connected  through  an  electric  clutch,  which  provides  for 
complete  decoupling  of  the  air  motor  when  free  rotation  of  the  sting  is  re¬ 
quired.  The  clutch  and  air  motor  can  also  be  used  for  braking  and  stopping 
the  sting  rotation. 

The  model  is  attached  to  the  sting  through  a  second  bearing  assembly 
which  allows  the  model  to  rotate  freely  with  respect  to  the  sting  (Figure  38). 
Alternately,  provision  is  made  for  locking  the  model  at  fixed  orientations 
with  respect  to  the  sting.  When  locked  to  the  sting  the  model  undergoes 
forced  lunar  motion. 

The  test  system  incorporates  a  series  of  stings  to  permit  testing  at 
different  angles  of  attack.  At  present,  15-,  30-,  and  45-degro**  stings  are 
available.  For  calibration  purposes,  a  special  zero-offset  sting  is  provided. 
The  zero-offset  sting  contributes  negligible  aerodynamic  damping,  but  does 
not  allow  the  model  to  roll  about  its  own  axis.  A  precision  angular  rate 
recording  device  is  provided  to  measure  the  sting  rotational  rate.  Model 
roll  behavior  must  be  recorded  optically  at  present, 

A  photograph  of  the  test  system  installed  in  the  Air  Force  Armament 
Laboratory  subsonic  wind  tunnel  is  shown  in  Figure  3‘).  In  the  photograph 
the  special  calibration  sting  is  installed. 


n 


Figure  37.  Dynamic -Wind -Tunnel -Model  Support  System 
Design  of  Rotating  Sting 


Figure  39.  Dynamic  Wind  Tunnel  Model  Support  System 
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R.  SYSTEM  OPERATION 

The  test  system  has  two  basic  modes  of  operation.  In  the  first  mode, 
the  model  is  driven  by  the  air  motor  to  a  pre-determined  coning  rate,  and 
then  disengaged  such  that  the  subsequent  motion  can  be  recorded.  Normally, 
this  will  be  a  deceleration  mode.  The  second  mode  of  operation  does  not 
require  the  drive  motor,  and  the  coning  is  initiated  by  the  inherent  autorota- 
tional  tendencies  of  the  model.  The  autorotation  can  be  produced  by  the 
magnus  moment  resulting  from  model  axial  spin  (model  spin  achieved  by  fin 
cant)  or  by  the  aerodynamic  side  momeht.  Aerodynamic  side  moments  can 
be  produced  by  asymmetric  roll  orientation  with  respect  to  the  angle  of 
attack  plane,  or  by  other  aerodynamic  asymmetries. 

To  preclude  autorotational  moments  from  model-sting  misalignment, 
the  boring  head  is  adjusted  so  that  the  sting  rotational  axis  and  the  model 
axis  of  symmetry  intersect  to  the  required  accuracy. 


C.  SYSTEM  UTILIZATION 

The  primary  purpose  of  the  test  system  is  the  measurement  of  the 
model  damping  for  circular-type  motions.  In  accordance  with  the  theory  of 
Reference  28,  it  can  be  anticipated  that  the  damping  coefficient  for  coning 
motion  will  differ  from  that  for  planar  pitching  motion  at  finite  angles  of 
attack.  To  date,  few  definitive  experiments  have  been  conducted  to  deter¬ 
mine  the  difference  between  the  circular  and  planar  motion  damping  coeffi¬ 
cients,  and  no  data  have  been  published  for  subsonic  conditions.  The 
circular  motion  damping  is  of  particular  interest  with  respect  to  the  S-Curve 
type  bomblet,  since  this  configuration  experiences  almost  pure  coning  motion 
in  free-flight. 

The  test  system  can  be  utilized  to  measure  magnus  and  side  mo¬ 
ments,  both  from  the  transient  and  steady- state  autorotational  coning 
motions. 

Roth  static  and  dynamic  roll  lock-in  tests  can  be  accomplished. 

Static  roll  lock-in  tests  are  accomplished  with  model  locked  to  the  sting  at 
specific  roll  orientations,  and  at  those  orientations  the  coning  and  side 
moment  characteristics  are  determined.  Dynamic  roll  lock-in  tests  are 
accomplished  for  those  models  which  have  inherent  aerodynamic  roll  lock-in 
tint*  to  body-fin  interference  effects.  The  tendency  for  the  lock-in  to  be 
stable  under  a  wide  range  of  coning  rates  and  fin  cant  angles  can  be  investi¬ 
gated. 


'['lie  test  system  also  provides  an  expedient  means  for  determining 
it  ',-i«d  roll  speed-up  and  roll  slow-down  at  large  angles  of  attack  and  the 
«  ffe<  t  «ii  coning  moti his  type  of  dynamic  behavior. 
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Procedure  for  Determination  of  Model  Circular  Motion  Damping 
The  circular  motion  damping  is  determined  from  time  histories  of  the  coning 
rate  decay,  in  a  manner  analagous  to  the  determination  of  the  roll  damping 
coefficient  in  single  degree-of-freedom  spin  decay  experiments.  Two  separ¬ 
ate  tares  may  be  required;  one  for  the  sting  bearing  friction  torque,  and  a 
second  for  the  offset  sting.  This  latter  tare  is  determined  by  the  use  of  the 
calibration  sting.  The  second  tare  is  not  required  if  the  model  is  to  always 
be  in  lunar  circular  motion,  but  in  general,  the  circular  motion  damping 
should  be  determined  for  various  model  roll  rates,  which  requires  the  off¬ 
set  sting. 

Using  the  nomenclature  shown  in  the  test  schematic.  Figure  40,  the 
equations  of  motion  for  each  model-sting  assembly  can  be  expressed  as: 


Ixfi  -  Tu  Si  -  Tq  =  0  (14) 

where 

Til  ’  *  <TU)„ 

To  =  (T„la  4  (T„)v 

and  subscript  a)  denotes  aerodynamic  torque  on  model 
v)  denotes  viscous  (or  bearing)  torque. 


The  quantity  Ix  is  the  total  moment  of  inertia  of  the  model  and  rotating 
sting  assembly,  taken  about  the  axis  of  sting  rotation. 


The  solution  of  equation  (14),  assuming  constantant  coefficient  is: 


+ 


)  exp  {Tw 


( t  -  t0  )/lx  I 


(IS) 


To  obtain  the  aerodynamic  torques,  the  system  is  first  operated 
with  the  model  and  offset- sting  section  removed  such  that  (T  ^)v  and 
(T0)v  can  be  evaluated.  Then  the  complete  model-sting  system  is  operated 
and  (Tq)  and  (T0)  are  obtained,  from  which  (Tq  )a  and  (T0)a  are 
readily  determined.  To  obtain  Tjj  and  T0  from  equation  ( 1 5),  a  lea 
squares  differential -correction  fitting  technique  is  employed,  as  desc: 
in  Reference  29.  The  possible  nonlinear  dependence  of  Tq  upon  fl  ono  ■*. 
ascertained  by  piecewise  fitting  the  .?  history. 


When  the  model  is  operated  with  the  offset  sting,  the  torque  ■  <  ''  .  ^ 
and  (Tu  )a  include  a  small  but  significant  contribution  from  the  stir.  i 
replacing  the  offset  sting  with  the  calibration  sting  (which  has  «tegltg>!>: 
aerodynamic  torque)  the  aerodynamic  torque  contribution  of  the  of.  -a 
sting  is  readily  determined.  Subtracting  the  offset  sting  torque  from  (T^)a 
and  (Tfl  )a  we  obtain  (T,  )a#  MODF1.  and  (Tt>  )a§  MODE1,. 


The  aerodynamic  torques  (Tfi  )a, MODEL  an<i  (Tc  )a> MODEL  are 
equated  to  the  model  spin  damping  and  yaw  damping  coefficients  by  the  fol¬ 
lowing  relationships: 

for  lunar  motion  (model  not  free  to  roll) 


^a,  MODEL 
1/2  p  Sd 


Cnr  ft  sin2  a  +  !1  cos  a  (16a) 


<To>a,  MODEL 
1/2  P  V2  Sd 


Cn  sin  a  *  Cc  cos  a 
no 


(16b) 


for  model  free  to  roll 


(Tft  ^a,  MODEL 
1/2  0  V2  Sd 


ft  sin2  a 


(16c) 


*to^,MODEL 
1/2P  V2  Sd 


( I6d) 


The  coefficients  Cj>  and  CnQ  arc  due  only  to  asymmetries  in  the 
model  and  flow  and  should  be  quite  small.  The  Coefficient  C  ^  will  usually 
be  evaluated  from  separate  tests  and  is  generally  small  compared  to  Cnr. 
Present  estimates  at  a  ’  1 5  degrees  give  CRr  =  -25.0  and  C  ^  =  -0.4. 

By  use  of  equations  (16a)  or  (16c)  Cn{,  can  be  evaluated  for  various 

values  of  angle  of  attack  and  ••  ,  These  data  can  then  be  used  to  determine 
possible  uonlinearities  with  respect  to  either  variable. 
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APPENDIX  I 

THEORETICAL  CONING  MOTION  ANALYSIS 
WITH  ROLL  DEPENDENT  SIDE  MOMENT 


The  fundamental  equations  describing  the  coning  motion  of  S-Curve 
type  bomblets  are  developed  in  Appendix  II  of  Reference  1.  Considered  here 
is  the  effect  of  a  roll  dependent  side  moment  on  the  stability  of  a  zero  coning 
solution. 

The  basic  equations  of  coning  motion  for  the  S-Curve  bomblet,  as 
originally  written,  do  not  include  the  effect  of  the  aerodynamic  roll  angle, 
i  .  Furthermore,  the  non-rotating  axes  system  which  is  utilized  does  not 
readily  permit  the  inclusion  of  body-fixed  aerodynamic  moments.  However, 
for  the  special  case  of  zero  spin,  the  effect  of  the  rotation  of  the  angle  of 
attack  plane  (variable  ^  )  and  the  aerodynamic  roll  angle  (  $  )  are  related 
simply  by 


U-l) 


The  special  coning  solution  may  now  be  written,  using  the  notation  of  Refer¬ 
ence  1.  as 


Sy  0 


f  ;  g  U  T  sin  f  1  CN  y  S  V  C„r  P  Sd*  V  ^  1  . 

-  ~  -  —  ♦  4  * 


g  pSd^  eos  i  sin  i-  (  *)  tfV^Sd 


»2  t? 


(1-2) 


inspection  of  the  order  «f  magnitude  of  the  last  two  terms  reveals 
that  loi*  a  six-fin  bomblet  a  steady •  state  solution  of  equation  (I*  2)  with 

0  requires  that  i  *  0,  “/h.  *  H,  ....  such  that  (tl  is  small, 
igure  l-l  illustrates  a  typical  variation  of  with  $  tot  the  basic 

S-Curve  configurations.  If  the  *ero  toning  solution  correspond*  to  J  »  0 
and  CRr  *  0*  then  (  $)  must  be  positive  and  vice  versa. 


To  determine  which,  if  any,  t»l  the  steady- state  *ero  coning  solutions 
arc  stable  solutions,  the  nature  of  the  characteristics  differential  equation 
for  the  small  perturbation  /*;•.; s  examined..  Th is  requires  that  the  term 
containing  M  )  be  made  a  lunation  of  .  This  can  be  accomplished 

by  considering  the  slope  of  versus  i  .  Thus,  equation  (1-2)  can  be 

f»- -w  ritten 


wemao  H3M  Kjuum  num. 


C$MsCSMt  sin  6$ 


1  aT  At  -  {  }  At 


Cnr  R  pSd  cos  Y 


(sin  $0  1  ) 


(1-3) 


d  CSM  l  2 

f  IT  3  PV  M  A* 


For  a  stable  solution,  one  of  the  requirements  is  that  the  coefficient  of  At 
be  positive.  Again,  inspecting  the  order  of  magnitude  of  the  two  terms 
which  represent  the  coefficients  of  At  .  it  is  found  that  the  d  CsM^  ^ 
term  is  much  larger  than  the  term  containing  Cnr.  Therefore,  a  stable 
motion  occurs  only  where  has  a  positive  slope.  Referring  to 

the  plot  of  versus  $  it  is  clear  that  i>  -  0  is  an  unstable  condition, 

and  that  the  bomblet  must  roll  approximately  $  =  1  30  degrees  to  where 
CSM  i®  small  and  has  a  positive  slope. 

The  coroila-y  to  the  above  is  that  a  change  in  sign  of  the  side  moment 
coefficient,  ,  should  make  the  bomblet  motion  stabie  at  <p  =  $  =  C. 

This  stability  hypothesis  has  been  confirmed  by  actual  6 -OOF  motion  sim¬ 
ulations  with  the  sign  of  the  side  moment  coefficient,  (^SM  i  *  purposely 
changed. 
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APPENDIX  II 

THEORETICAL  PREDICTION  OF  CLUSTER  SPREAD  VELOCITY 


A  tV>f»n r-ettcal  estimate  of  the  maximum  cluster  spread  velocity  is 
made,  by  considering  the  motion  of  a  single  bomblet  initially  at  the  leading 
edge  of  the  cluster.  The  forces  acting  on  the  bomblet  are  computed  from 
the  potential  flow  about  a  solid  body  assumed  to  represent  the  core  of  the 
cluster.  This  model  is  predicated  upon  observation  of  actual  bomblet 
clusters,  where  it  appe?.rs  that  the  spreading  is  taking  place  primarily  at 
the  leading  edge  of  the  cluster. 

As  a  first  approach  at  a  solution,  consider  a  core  of  spherical  shape 
with  a  diameter  closely  related  to  that  of  the  dispenser.  The  physical  model 
of  the  problem  is  like  that  sketched  below. 


The  aerodynamic  forces  acting  me.  the  bomblet  *  t  *  aU>s  4  the  indicated 
lateral  trajectory  can  now  be  evaluated,  starting  from  an  tn-tia*  point 
slightly  offset  from  the  dividing  streamline  and  extending  Utefaliv  until 
near  uniform  flow  is  encountered.  The  selection  of  this  parti*  -la;  trajec¬ 
tory  is  equivalent  to  the  assumption  that  the  bomblet  and  the  cluster  have 
the  same  decele-ation. 

The  computation  of  the  flow  direction  and  magnitude  at  each  point 
along  this  trajectory  is  easily  ci-mplishtd  using  a  known  polemic!  flow 
solution,  (for  example,  pp  ?*»,  "Aerodynamics  of  Wings  and  J  odies. 

Ashley  and  Landahl.  Addison-Wesley  Press). 
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Since  the  direction  and  magnitude  of  the  flow  at  each  point  along  the 
radial  trajectory  is  known,  the  aerodynamic  lateral  force  can  also  be  com¬ 
puted  at  each  point.  To  begin,  consider  only  the  bomblet  drag  force,  and 
its  component  in  the  lateral  direction,  Fy.  The  force,  Fy,  when  integrated 
as  a  function  of  the  radial  deflection,  represents  the  work  done  on  the  bomb- 
let,  and  it  can  be  equated  to  the  radial  kinetic  energy  of  the  bomblet.  This 
in  turn  can  be  evaluated  in  terms  of  the  radial  velocity. 

Defining  the  parameters  in  accordance  with  the  sketch  below,  the 
radial  velocity  can  be  derived  as  follows: 


The  work  done  by  the  drag  force  can  be  expressed  as 

rv 

Work  /  ly  dY 


m-n 


which  can  be  re-written  as 


f: 


Work  -  /  {  ?U)  CD  S  V2  { ■—*)  d  Y 


Introducing  the  new  variable  Y/K,  and  relating  the  velocities  to 

.Y/R 


Work  i  4  C 

6 


*’su-R/  [ufM*] 


(U-2) 
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Using  the  potential  flow  solution  for  a  sphere  at  X/R  =  1.0,  the  quantity 
v  V 

U“7  has  been  evaluated  and  is  plotted  in  Figure  II -1  as  a  function  of  Y/R. 

00 

It  is  seen  that  most  of  the  work  done  on  the  bomblet  occurs  with  Y/R  <  3.  0. 
A  very  good  approximation  of  the  total  work  is  achieved  by  integrating  to 
Y/R  =5.0. 


By  equating  the  work  done  on  the  bomblet  to  its  radial  kinetic  energy, 
a  relationship  for  the  radial  velocity  is  obtained  in  the  form 


1  2  ? 

~  mn  Vo  (  p  / 2 )  CnSU  R 


2  mR  R 


f  Y/R 

vV  ] 

d  fl] 

/ 

U2 

CU  ^ 

d  I  R 

L  J 

Y/R 


or 


U 


IP  R 


1W/CD  S)B 


vV 


U 


_Y 

R 


1 

'2 


(II -2) 


The  integral  can  be  evaluated  from  Figure  II- 1  and  has  a  value  of  approxi 
mately  0.39.  Therefore  approximately 


■R 


U 


where 


1/2 


0.39 


K  P  R 


(\V/CnS) 


J 


( II  -4) 


R  Cluster  core  radius 

<  V\ /t’ljSijj  bomblet  ballistic  coefficient. 


Thus,  the  radial  velocity  imparted  to  a  bomblet  displaced  from  the  leading 
edge  of  the  cluster  is 


1)  Directly  proportional  to  the  cluster  (event)  velocity. 

2)  Proportional  to  the  square -rout  of  the  cluster  sut  . 

3'  Inversely  proportional  to  the  square- root  of  th»-  bomblet 
ballistic  coefficient. 

Equation  01  *4)  was  evaluated  for  both  .md  A  PAM 

clusters,  assuming  that  R  was  the  same  as  the  dispenser  radius.  For 
tue  CI3U -70/13  the  radial  velocity  is  determined  to  i>e 

V,.  0.018  U 

VM  =  0.013 


and  for  A  HAM 


*/R  *  1.0 


y'R 


Figure  1 1  - 1 .  lateral  Force  Parameter  for  Nomblet  Separating 
from  Cluster 
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Test  data  (Figure  33)  indicate  a  value  of  V^ /Uro  =  0.03  at  U^  =  900  feet/ 
second.  Thus  the  method  appears  to  give  a  correct  order  of  magnitude  pre¬ 
diction,  although  the  predicted  radial  velocity  is  somewhat  low.  This  may 
be  due  to  selecting  R  as  the  dispenser  radius;  possibly  the  cluster  can 
expand  to  larger  size  and  still  behave  as  a  quasi-solid  core. 

The  above  theory  confirms  the  findings  of  Reference  19  that  the 
spread  velocity  is  a  linear  function  of  the  event  velocity.  The  theoretical 
dependence  of  the  radial  velocity  on  the  bomblet  ballistic  coefficient  is  a 
new  result.  It  will  be  of  future  interest  to  see  i<‘  this  can  be  confirmed  by 
test  data. 
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Pom?"  carlo  analyses.  Improved  aerodynamic  data  packages  for  both  the  3-Curve 
-ci-ro  bomb  lets  are  provided,  and  results  of  additional  wind  tunnel 
tests  rtf-  discussed •  Monte  Carlo  simulations  show  that  for  realistic  flight  en¬ 
vironments.  both  tie-  3-Curve  and  rol 1 -through-zero  bomblets  can  achieve  large  and 
uniform  impact  pattern'...  For  a  representative  dispenser  HOB  of  2000  feet,  opening 
velocity  of  900  fpo ,  and  1*5  degrees  flight  path  angle,  impact  pattern  widths  of 
HuO  and  000  feet  are  computed  for  the  0-Curve  ami  roll-through-zero  type  bomblets , 
respectively.  The  effects  of  transonic  and  supersonic  delivery,  large  static  mass 
unbalance,  nose-roughness  asymmetry,  and  intentional  fin  cant  on  the  S-Curve 
bomblet  flight  characteristics  and  dispersion  are  investigated  in  detail.  Flight 
test  dispersion  data  for  the  S-Curve  bomblet ,  as  obtained  from  air-gun  launched 
models ,  is  compared  with  results  from  analytical  simulations.  A  dynamic  wind 
tnnncT  model  support  system,  designed  for  testing  the  S-Curve  bomblet,  is  de¬ 
scribe.!. 
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